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FUTURE  LAND  USE  SCENARIO  .  ADULTS  AND  CHILDREN 


TABLE  P-68  CHRONIC  HAZARDOUS  INDEX  ESTIMATE  -  SURFACE  WATER  INGESTION 
FUTURE  LAND  USE  SCENARIO  -  ADULTS  AND  CHILDREN 


Noncardnogehc  Effects 

Site  1 ,  Former  POL  Facility,  Alpena 
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TABLE  P-60  SUBCHRONIC  HAZARDOUS  INDEX  ESTIMATE  -  SOIL  INHALATION 
FUTURE  LAND  USE  SCENARIO  -  EXCAVATION  WORKER 


Groundwater  Solute  Transport  Model  Data  •  Site  1 
Theoretical  Background 

A  two-dimensional  Method  of  Characteristics  (MOO  solute  transport  model  (Konikow  and 
Bredehoeft,  1989)  was  used  for  preliminary  examination  of  contaminant  migration  within 
the  shallow  aquifer  beneath  the  Alpena  CRTC.  The  model  is  designed  to  calculate 
transient  changes  in  solute  concentrations  within  groundwater  by  simultaneously  solving 
partial  differential  equations  describing  groundwater  flow  and  transport  and  computes  the 
change  in  a  chemicals  concentration  over  time.  Changes  in  chemical  concentrations  over 
time  are  caused  by  the  processes  of  convective  transport,  hydrodynamic  dispersion,  and 
mixing  from  fluid  sources.  This  model  couples  the  groundwater  flow  equation  with  solute 
transport  equations. 

The  flow  equation  can  be  approximated  by  an  implicit  finite-difference  equation.  The 
model  area  is  discretized  into  a  rectangular  grid  with  each  square  being  a  node.  The  finite 
difference  equation  is  solved  numerically  for  each  node  in  the  grid  using  an  iterative 
alternating-direction  implicit  (ADI)  procedure. 

After  the  hydraulic  head  distribution  is  calculated,  the  velocity  of  groundwater  flow  can  be 
computed  at  each  node.  The  expression  for  average  velocity  of  groundwater  can  be 
derived  from  Darcy's  law.  The  groundwater  velocity  at  each  node  is  calculated  utilizing  an 
explicit  finite-difference  approximation  of  Darcy's  law.  The  computer  program  uses  an 
alternating-direction  implicit  procedure  to  solve  a  finite-difference  approximation  to  the 
groundwater  flow  equation,  and  it  uses  the  method  of  characteristics  (MOC)  to  solve  the 
solute  transport  equation.  MOC  uses  a  particle  tracking  procedure  to  represent  convective 
transport  and  a  two-step,  explicit  procedure  to  solve  a  finite-difference  equation  that 
describes  the  effects  of  hydrodynamic  dispersion,  fluid  sources  and  sinks,  and  divergence 
of  velocity. 

A  number  of  assumptions  are  inherent  in  the  solute  transport  model: 

1.  Darcy's  law  is  valid  and  hydraulic  head  gradients  are  the  only  significant  driving 
mechanism  for  fluid  flow. 

2.  The  porosity  and  hydraulic  conductivity  of  the  aquifer  are  constant  with  time,  and 
porosity  is  uniform  in  space. 

3.  Gradients  of  fluid  density,  viscosity,  and  temperature  do  not  affect  the  velocity 
distribution. 

4.  No  chemical  reactions  occur  that  affect  the  concentration  of  the  solute,  the  fluid 
properties,  or  the  aquifer-properties. 

5.  Ionic  and  molecular  diffusion  are  negligible  contributors  to  the  total  dispersive  flux. 

6.  Vertical  variations  in  head  and  concentrations  are  negligible. 
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7.  The  aquifer  is  homogeneous  and  isotropic  with  respect  to  the  coefficients  of 
longitudinal  and  transverse  dispersivity. 


Transport  Model  Input 

A  model  grid  of  32  columns  by  19  rows  with  a  250  foot  lateral  spacing  was  used. 
Specified  head  cells  were  used  at  nodes  corresponding  to  the  South  Branch  of  the  Thunder 
Bay  River,  at  nodes  along  the  eastern  boundary  of  the  model  grid  area,  and  also  at  the 
sinkhole  in  the  northeastern  portion  of  the  model.  Groundwater  elevations  measured 
during  September,  1993  were  used  as  initial  input  into  the  transport  model. 

Hydraulic  conductivity  values  were  calculated  from  slug  tests  performed  at  Alpena  CRTC 
(Engineering  Science,  1989;  Earth  Technology,  1994).  Values  of  hydraulic  conductivity 
range  from  12  feet/day  at  Site  4  to  278  feet/day  at  Site  3. 

Aquifer  thickness  values  were  obtained  from  drilling  records  of  monitoring  wells  and  soil 
borings  obtained  from  the  SI  and  Rl  field  activities.  Values  listed  are  from  logs  in  which 
the  thickness  of  the  shallow  aquifer  was  clearly  discernible,  and  ranged  from  20  feet  at 
Site  5  to  65  feet  at  Site  8.  Transmissivity  values  were  calculated  by  multiplying  the 
calculated  hydraulic  conductivity  values  by  the  aquifer  thickness.  Transmissivity  ranges 
from  420  ft^/day  at  TF4-MW3  to  15,290  ft^/day  at  CG3-MW5. 

Monitoring  of  the  discharge  of  springs  into  the  sinkhole  was  performed  during  the  SI 
(Engineering  Science,  1990)  and  an  estimate  of  approximately  18,000  gallons  of  water  per 
day  discharging  into  the  sinkhole  was  calculated.  In  order  to  obtain  a  numerical  estimate 
of  discharge  into  the  sinkhole  for  the  model,  MODFLOW  (McDonald  and  Harbaugh,  1988), 
a  3-  dimensional  finite-difference  groundwater  flow  model  was  used.  MODFLOW  was 
used  because  of  its  ability  to  simulate  the  effect  of  head-dependent  groundwater  flow  into 
a  groundwater  sink  (i.e.  the  sinkhole).  This  package  was  not  available  in  MOO.  The  same 
model  parameters  and  boundary  conditions  were  used  within  MODFLOW  as  in  MOC. 

Based  upon  hydraulic  head  data  collected  in  September  1993,  discharge  from  the  shallow 
aquifer  into  the  sinkhole  is  approximately  30,000  gallons  per  day. 

The  dispersivity  of  an  aquifer  in  two  dimensions  is  described  by  the  longitudinal  dispersion, 
the  transverse  dispersion  and  the  ratio  of  the  two  (Fetter,  1993).  As  a  contaminant  plume 
moves  further  from  its  initial  location  within  the  aquifer  by  advection  with  the  groundwater 
flow,  the  plume  spreads.  The  spreading  in  the  direction  of  groundwater  flow  is  the 
longitudinal  dispersion,  the  spreading  in  the  direction  perpendicular  to  the  groundwater 
flow  is  known  as  the  transverse  dispersion  (Fetter,  1993).  The  values  of  the  dispersion 
coefficients  are  typically  derived  via  bench  scale  tests,  aquifer  tests,  or  calibration  of 
contaminant  transport  models.  Since  no  data  presently  exists  describing  dispersivity 
within  the  shallow  aquifer  beneath  the  Alpena  CRTC  and  insufficient  data  exists  to  allow 
for  derivation  of  dispersivity  via  model  calibration,  moderate  values  of  100  feet  for 
longitudinal  dispersivity  and  30  feet  for  transverse  dispersivity  were  chosen  (Gillham  and 
Cherry,  1982).  A  more  complete  description  of  the  model  is  given  in  the  report. 

Preliminary  Groundwater  Modeling  Effort,  Earth  Technology,  August  1993. 
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Model  Calibration: 


The  groundwater  flow  model  was  calibrated  with  respect  to  the  September  1993 
groundwater  elevation  measurements.  Calibration  of  the  groundwater  flow  model  was 
accomplished  by  defining  a  set  of  parameters,  boundary  conditions,  and  stresses  that 
produce  simulated  heads  and  fluxes  that  match  field-measured  values  within  a 
preestablished  range  of  error  (Anderson  and  Woessner,  1992).  In  order  to  match  field 
measured  values  for  hydraulic  head  as  determined  during  September  1 993,  a  few 
modifications  were  made  to  the  preexisting  groundwater  flow  model.  These  changes 
included  updating  the  initial  head  array,  modeling  the  sinkhole  as  a  constant  head  cell  to 
account  for  the  large  gradient  changes  in  the  vicinity  of  the  sinkhole  and  including  recharge 
to  the  model  at  a  rate  of  9  inches  per  year  over  the  whole  model  area.  By  adjusting  these 
parameters,  an  acceptable  level  of  calibration  was  achieved.  An  acceptable  level  of 
calibration  was  defined  as  a  root  mean  squared  error  (RMS)  of  less  than  2  feet.  The  RMS, 
or  the  standard  deviation  is  the  average  squared  difference  in  measured  and  simulated 
heads  and  is  given  by  the  equation: 


n 


RMS=[l/nY^ 


o.s 


n  =  number  of  wells 
h„  =  measured  head 
h,  =  model  simulated  head 


Certain  portions  of  the  model  may  have  values  above  the  goal  of  2  feet  while  others  fall 
much  below  this  value.  The  RMS  represents  the  average  error  present  in  the  model. 
The  following  provides  a  summary  of  the  final  calibrated  heads  for  the  flow  model. 


Well  # 

Row,  Col 

Measured 

Head  (ft) 

Simulated 

Head  (ft) 

Measured  - 
Simulated 

Head  (ft) 

MP2MW1 

10,25 

679.69 

676.68 

3.01 

MP2MW3 

10,22 

675.51 

675.46 

0.05 

MP2MW4_5* 

11,23 

675.34 

675.70 

-0.36 

MP2MW6 

11,20 

674.86 

674.52 

0.34 

CG3MW1 

6,24 

677.38 

676.69 

0.69 

CG3MW2 

7,22  - 

676.29 

675.78 

0.51 

CG3MW3 

9,23 

676.50 

675.98 

0.52 

CG3MW4_5* 

8,23 

676.41 

676.08 

0.33 
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Well  # 

Row,  Col 

Measured 

Head  (ft) 

Simulated 

Head  (ft) 

Measured  • 
Simulated 

Head  (ft) 

CG3MW7 

8,20 

675.64 

674.80 

0.84 

TF4MW1 

7,9 

667.23 

668.83 

-1.60 

TF4MW2 

8,10 

659.61 

664.22 

-4.61 

TF4MW3_4* 

9,10 

658.21 

660.06 

-1.85 

SF5MW1 

12,6 

674.15 

671.34 

2.81 

SF5MW2 

14,5 

675.32 

672.67 

2.65 

SF5MW3_4* 

13,6 

676.46 

671.64 

4.82 

SF5MW6 

13,5 

674.26 

672.40 

1.86 

LF6MW1 

14,9 

672.1 

670.40 

1.70 

LF6MW2 

14,8 

672.68 

670.72 

1.96 

LF6MW3 

14,10 

671.93 

670.21 

1.72 

LF6MW4 

16,7 

672.75 

671.91 

0.84 

LF6MW5 

15,7 

673.07 

671.56 

1.51 

LF6MW6 

13,10 

671.17 

669.67 

1.50 

LF6MW8 

15,9 

673.12 

670.78 

2.34 

HN8MW1 

5,22 

676.96 

675.93 

1.03 

HN8MW2 

6,19 

675.31 

674.50 

0.81 

HN8MW3_4* 

7,21 

676.01 

675.35 

0.66 

RTS  MW  1 

6,16 

673.06 

672.78 

0.28 

RT9MW2 

7,14 

668.21 

670.81 

-2.60 

RT9MW3 

9,15 

670.72 

671.26 

-0.54 

RT9MW4_5* 

8,14 

667.47 

670.32 

-2.85 

RT9MW6 

8,16 

670.58 

672.33 

-1.75 

S1MW2 

13,26 

677.39 

676.63 

0.76 

S1MW3 

13,27 

677.15 

676.98 

0.17 

S1MW1 1 

15,24 

675.72 

675.67 

0.05 

S1MW12 

16,25 

674.55 

675.94 

-1.39 
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Well  # 

Row,  Col 

Measured 

Head  (ft) 

Simulated 

Head  (ft) 

Measured  • 
Simulated 

Head  (ft) 

S1MW13 

14,24 

675.21 

675.76 

-0.55 

S1MW14 

14,25 

673.92 

676.14 

-2.22 

MP2MW2 

12,24 

675.57 

675.96 

-0.39 

*  Indicates  that  more  than  one  well  is  present  in  each  node  and  an  average  value  for 
hydraulic  head  was  used. 

Sum  of  Squared  Residuals  =  128.78/38  =  3.3891 
Root  Mean  Squared  Error  =  1 .84 


It  should  also  be  noted  that  the  model  was  calibrated  with  respect  to  the  September  1993 
water  level  data  and  should  only  be  considered  calibrated  with  respect  to  this  data.  More 
information  on  the  water  level  fluctuation  through  time  would  be  needed  to  perform  a 
transient  calibration.  The  model  was  not  calibrated  with  respect  to  concentration  data,  but 
only  with  respect  to  hydraulic  head. 

Model  Assumptions  and  Limitations 

*  The  model  domain  consisted  only  of  the  shallow  unconfined  aquifer  (i.e.  one  layer). 

*  The  initial  head  data  input  to  the  transport  model  are  results  of  measurements  taken  in 
September  1993. 

*  Initial  concentrations  of  compounds  are  results  of  the  Round  IV  sampling  event  which 
was  conducted  from  July  to  September  1993. 

*  Hydraulic  conductivity  values  are  the  result  of  slug  tests  performed  in  November,  1987 
and  September  1993. 

*  The  model  was  calibrated  with  respect  to  hydraulic  head  using  September  1993  water 
level  data  and  should  only  be  considered  calibrated  with  respect  to  September  1993 
water  level  data. 

*  The  flow  model  was  assumed  to  be  at  steady-state  with  respect  to  hydraulic  head. 


site  1 


The  following  compounds  were  present  at  Site  1:  Benzene,  Bromodichloromethane, 
Dibromochloromethane,  1 ,4-  Dichlorobenzene,  Styrene,  Antimony,  Carbon  Tetrachloride, 
and  Dibenzofuran.  Concentrations  input  to  the  model  are  the  result  of  groundwater 
sampling  conducted  during  the  July  to  September  1993  Rl  field  events. 

For  Site  1 ,  several  scenarios  were  modeled.  These  included  monitoring  chemical 
concentrations  through  time  along  the  South  Branch  of  the  Thunder  Bay  River  and  at 
Production  Well  #3  for  chemicals  observed  in  several  wells  at  Site  1.  Wells  SI  MW  14, 
S1MW1 1,  S1MW6,  and  Production  Well  #3  all  had  chemicals  present  above  MDNR  Type 
A  or  B  cleanup  criteria. 

The  following  scenarios  were  modeled  for  S1MW14:  Benzene  at  1 1.5  ug/1, 
Bromodichloromethane  at  0.78  ug/l,  Dibromochlormethane  at  2.1  ug/l,  1,4  - 
Dichlorobenzene  at  0.08  ug/l,  and  Styrene  at  a  concentration  of  0.13  ug/I.  Concentrations 
were  input  to  the  node  corresponding  to  SI  MW  14  and  concentrations  through  time  were 
monitored  along  the  South  Branch  of  the  Thunder  Bay  River. 

Wells  S1MW1 1  and  S1MW6  contained  elevated  levels  of  Antimony  at  a  concentration  of 
39.2  ug/l.  This  concentration  was  input  to  the  model  and  the  corresponding 
concentrations  of  Antimony  through  time  were  monitored  along  the  South  Branch  of  the 
Thunder  Bay  River. 

Additionally,  Carbon  Tetrachloride  and  Dibenzofuran  were  present  in  Production  Well  ttZ  at 
concentrations  of  1.2  ug/l  and  2.5  ug/l,  respectively.  These  values  were  input  to  the 
model  and  the  concentrations  were  monitored  through  time  along  the  South  Branch  of  the 
Thunder  Bay  River. 
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Appendix  Q:  Site  2  Risk  Assessment 
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Table  Q-2A 

Data  Utilized  to  Calculate  Reasonable  Maximum  Exposure  Concentrations  for  Soil  -  Site  2 

MIANG,  Alpena  CRTC,  Alpena,  Ml 
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TABLE  Q-2  DAILY  INTAKE  -  ADULT 

INGESTION  OF  SHALLOW  GROUNDWATER 
FUTURE  LAND  USE  SCENARIO 
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TABLE  Q-6  DAILY  INTAKE  -  ADULT 


TABLE  Q-8  DAILY  INTAKE  -  ADULT 
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ASSESSMENT  PARAMETERS  -  INHALATION  OF  GROUNDWATER 
I,  Alpena  CRTC,  Alpena,  Ml 


TABLE  Q-1 3  DAILY  INTAKE- ADULT 

INHALATION  OF  SHALLOW  GROUNDWATER 
FUTURE  LAND  USE  SCENARIO 


TABLE  Q-15  DAILY  INTAKE  -  ADULT 

INHALATION  OF  GROUNDWATER  DEEP  AQUIFER  PRODUCTION  WELLS 
FUTURE  LAND  USE  SCENARIO 
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Carcinogenic  Effect  Exposure  Assessment 
Site  2,  Motor  Pool,  Alpena  CRTC,  Alpena  Ml 
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TABLE  Q-19  GROUNDWATER  INHALATION  MODEL  CALCULATIONS 
Site  2,  Motor  Pool,  Alpena  CRTC,  Alpena,  Ml 


TABLE  Q-21  DAILY  INTAKE  -  ADULT 

SHALLOW  GROUNDWATER  DERMAL  CONTACT 
CURRENT  LAND  USE  SCENARIO 


Q-15 


TABLE  Q-23  DAILY  INTAKE  -  ADULT 

SHALLOW  GROUNDWATER  DERMAL  CONTACT 
CURRENT  LAND  USE  SCENARIO 


TABLE  Q-27  DAILY  INTAKE  -  ADULT 

GROUNDWATER  DERMAL  CONTACT  DEEP  AQUIFER  PRODUCTION  WELLS 
FUTURE  LAND  USE  SCENARIO 


TABLE  Q^O  CHRONIC  HAZARDOUS  INDEX  ESTIMATE -SHALLOW  GROUNDWATER  INGESTION 
FUTURE  LAND  USE  SCENARIO  -  ADULTS  AND  CHILDREN 


TABLE  Q^1  CANCER  ESTIMATE  -  GROUNDWATER  INGESTION  DEEP  AQUIFER  PRODUCTION  WELLS 
FUTURE  LAND  USE  SCENARIO  -  ADULTS  AND  CHILDREN 


TABLE  Q-32  CHRONIC  HAZARDOUS  INDEX  ESTIMATE  -  GROUNDWATER  INGESTION  DEEP  AQUIFER  PRODUCTION  WELLS 
FUTURE  LAND  USE  SCENARIO  -  ADULTS  AND  CHILDREN 


TABLE  Q-33  CANCER  ESTIMATE  -  SHALLOW  GROUNDWATER  INHALATION 
FUTURE  LAND  USE  SCENARIO  -  ADULTS  AND  CHILDREN 


TABLE  Q-34  CHRONtC  HAZARDOUS  INDEX  E3TWATE  -  SHALLOW  GROUNDWATER  INHALATION 
FUTURE  LAND  USE  SCENARIO  -  ADULTS  AND  CHILDREN 


TABLE  Q-38  CHRONIC  HAZARDOUS  INDEX  ESTWIATE  •  GROUNDWATER  INHALATION  DEEP  AQUIFER  PRODUCTION  WELLS 
FUTURE  LAND  USE  SCENARIO  -  ADULTS  AND  CHILDREN 


TABLE  Q-37  CANCER  ESTIMATE  -  DERMAL  CONTACT  WITH  SHALLOW  GROUNDWATER 
FUTURE  LAND  USE  SCENARIO  -  ADULTS  AND  CHILDREN 


TABLE  Q-38  CHRONIC  HAZARDOUS  INDEX  ESTIMATE  -  DERMAL  CONTACT  WITH  SHALLOW  GROUNDWATER 
FUTURE  LAND  USE  SCENARIO  -  ADULTS  AND  CHILDREN 


TABLE  Q-39  CANCER  ESTIMATE  •  DERMAL  CONTACT  WITH  GROUNDWATER  DEEP  AQUIFER  PRODUCTION  WELLS 
FUTURE  LAND  USE  SCENARIO  -  ADULTS  AND  CHILDREN 


TABLE  Q^O  CHRONIC  HAZARDOUS  INDEX  ESTIMATE  -  DERMAL  CONTACT  WITH  GROUNDWATER  DEEP  AQUIFER  PRODUCTION  WELLS 
FUTURE  LAND  USE  SCENARIO  -  ADULTS  AND  CHILDREN 


Groundwater  Solute  Transport  Model  Data  -  Site  2 
Theoretical  Background 

A  two-dimensional  Method  of  Characteristics  (MOC)  solute  transport  model  (Konikow  and 
Bredehoeft,  1989)  was  used  for  preliminary  examination  of  contaminant  migration  within 
the  shallow  aquifer  beneath  the  Alpena  CRTC.  The  model  is  designed  to  calculate 
transient  changes  in  solute  concentrations  within  groundwater  by  simultaneously  solving 
partial  differential  equations  describing  groundwater  flow  and  transport  and  computes  the 
change  in  a  chemicals  concentration  over  time.  Changes  in  chemical  concentrations  over 
time  are  caused  by  the  processes  of  convective  transport,  hydrodynamic  dispersion,  and 
mixing  from  fluid  sources.  This  model  couples  the  groundwater  flow  equation  with  solute 
transport  equations. 

The  flow  equation  can  be  approximated  by  an  implicit  finite-difference  equation.  The 
model  area  is  discretized  into  a  rectangular  grid  with  each  square  being  a  node.  The  finite 
difference  equation  is  solved  numerically  for  each  node  in  the  grid  using  an  iterative 
alternating-direction  implicit  (ADI)  procedure. 

After  the  hydraulic  head  distribution  is  calculated,  the  velocity  of  groundwater  flow  can  be 
computed  at  each  node.  The  expression  for  average  velocity  of  groundwater  can  be 
derived  from  Darcy's  law.  The  groundwater  velocity  at  each  node  is  calculated  utilizing  an 
explicit  finite-difference  approximation  of  Darcy's  law.  The  computer  program  uses  an 
alternating-direction  implicit  procedure  to  solve  a  finite-difference  approximation  to  the 
groundwater  flow  equation,  and  it  uses  the  method  of  characteristics  (MOC)  to  solve  the 
solute  transport  equation.  MOC  uses  a  particle  tracking  procedure  to  represent  convective 
transport  and  a  two-step,  explicit  procedure  to  solve  a  finite-difference  equation  that 
describes  the  effects  of  hydrodynamic  dispersion,  fluid  sources  and  sinks,  and  divergence 
of  velocity, 

A  number  of  assumptions  are  inherent  in  the  solute  transport  model. 

1.  Darcy's  law  is  valid  and  hydraulic  head  gradients  are  the  only  significant  driving 
mechanism  for  fluid  flow. 

2.  The  porosity  and  hydraulic  conductivity  of  the  aquifer  are  constant  with  time,  and 
porosity  is  uniform  in  space. 

3.  Gradients  of  fluid  density,  viscosity,  and  temperature  do  not  affect  the  velocity 
distribution. 

4.  No  chemical  reactions  occur  that  affect  the  concentration  of  the  solute,  the  fluid 
properties,  or  the  aquifer'properties. 

5.  Ionic  and  molecular  diffusion  are  negligible  contributors  to  the  total  dispersive  flux. 
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6.  Vertical  variations  inTiead  and  concentrations  are  negligible. 

7.  The  aquifer  is  homogeneous  and  isotropic  with  respect  to  the  coefficients  of 
longitudinal  and  transverse  dispersivity. 


Transport  Model  Input 

A  model  grid  of  32  columns  by  1 9  rows  with  a  250  foot  lateral  spacing  was  used. 

Specified  head  cells  were  used  at  nodes  corresponding  to  the  South  Branch  of  the  Thunder 
Bay  River,  at  nodes  along  the  eastern  boundary  of  the  model  grid  area,  and  also  at  the 
sinkhole  in  the  northeastern  portion  of  the  model.  Groundwater  elevations  measured 
during  September,  1993  were  used  as  initial  input  into  the  transport  model. 

Hydraulic  conductivity  values  were  calculated  from  slug  tests  performed  at  Alpena  CRTC 
(Engineering  Science,  1989;  Earth  Technology,  1994).  Values  of  hydraulic  conductivity 
range  from  1 2  feet/day  at  Site  4  to  278  feet/day  at  Site  3. 

Aquifer  thickness  values  were  obtained  from  drilling  records  of  monitoring  wells  and  soil 
borings  obtained  from  the  SI  and  Rl  field  activities.  Values  listed  are  from  logs  in  which 
the  thickness  of  the  shallow  aquifer  was  clearly  discernible,  and  ranged  from  20  feet  at 
Site  5  to  65  feet  at  Site  8.  Transmissivity  values  were  calculated  by  multiplying  the 
calculated  hydraulic  conductivity  values  by  the  aquifer  thickness.  Transmissivity  ranges 
from  420  ft^/day  at  TF4-MW3  to  1 5,290  ft^/day  at  CG3-MW5, 

Monitoring  of  the  discharge  of  springs  into  the  sinkhole  was  performed  during  the  SI 
(Engineering  Science,  1990)  and  an  estimate  of  approximately  18,000  gallons  of  water  per 
day  discharging  into  the  sinkhole  was  calculated.  In  order  to  obtain  a  numerical  estimate 
of  discharge  into  the  sinkhole  for  the  model,  MODFLOW  (McDonald  and  Harbaugh,  1988), 
a  3-  dimensional  finite-difference  groundwater  flow  model  was  used.  MODFLOW  was 
used  because  of  its  ability  to  simulate  the  effect  of  head-dependent  groundwater  flow  into 
a  groundwater  sink  (i.e.  the  sinkhole).  This  package  was  not  available  in  MOC.  The  same 
model  parameters  and  boundary  conditions  were  used  within  MODFLOW  as  in  MOC. 

Based  upon  hydraulic  head  data  collected  in  September  1993,  discharge  from  the  shallow 
aquifer  into  the  sinkhole  is  approximately  30,000  gallons  per  day. 

The  dispersivity  of  an  aquifer  in  two  dimensions  is  described  by  the  longitudinal  dispersion, 
the  transverse  dispersion  and  the  ratio  of  the  two  (Fetter,  1993).  As  a  contaminant  plume 
moves  further  from  its  initial  location  within  the  aquifer  by  advection  with  the  groundwater 
flow,  the  plume  spreads.  The  spreading  in  the  direction  of  groundwater  flow  is  the 
longitudinal  dispersion,  the  spreading  in  the  direction  perpendicular  to  the  groundwater 
flow  is  known  as  the  transverse  dispersion  (Fetter,  1993).  The  values  of  the  dispersion 
coefficients  are  typically  derived  via  bench  scale  tests,  aquifer  tests,  or  calibration  of 
contaminant  transport  models.  Since  no  data  presently  exists  describing  dispersivity 
within  the  shallow  aquifer  beneath  the  Alpena  CRTC  and  insufficient  data  exists  to  allow 
for  derivation  of  dispersivity  via  model  calibration,  moderate  values  of  100  feet  for 
longitudinal  dispersivity  and  30  feet  for  transverse  dispersivity  were  chosen  (Gillham  and 
Cherry,  1982).  A  more  complete  description  of  the  model  is  given  in  the  report. 
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Preliminary  Groundwater  Modeling  Effort,  Earth  Technology,  August  1993. 


Model  Calibration: 

The  groundwater  flow  model  was  calibrated  with  respect  to  the  September  1993 
groundwater  elevation  measurements.  Calibration  of  the  groundwater  flow  model  was 
accomplished  by  defining  a  set  of  parameters,  boundary  conditions,  and  stresses  that 
produce  simulated  heads  and  fluxes  that  match  field-measured  values  within  a 
preestablished  range  of  error  (Anderson  and  Woessner,  1992).  In  order  to  match  field 
measured  values  for  hydraulic  head  as  determined  during  September  1993,  a  few 
modifications  were  made  to  the  preexisting  groundwater  flow  model.  These  changes 
included  updating  the  initial  head  array,  modeling  the  sinkhole  as  a  constant  head  cell  to 
account  for  the  large  gradient  changes  in  the  vicinity  of  the  sinkhole  and  including  recharge 
to  the  model  at  a  rate  of  9  inches  per  year  over  the  whole  model  area.  By  adjusting  these 
parameters,  an  acceptable  level  of  calibration  was  achieved.  An  acceptable  level  of 
calibration  was  defined  as  a  root  mean  squared  error  (RMS)  of  less  than  2  feet.  The  RMS, 
or  the  standard  deviation  is  the  average  squared  difference  in  measured  and  simulated 
heads  and  is  given  by  the  equation: 


i-1 


n  =  number  of  wells 
hn,  =  measured  head 
h,  =  model  simulated  head 


Certain  portions  of  the  model  may  have  values  above  the  goal  of  2  feet  while  others  fall 
much  below  this  value.  The  RMS  represents  the  average  error  present  in  the  model. 

The  following  provides  a  summary  of  the  final  calibrated  heads  for  the  flow  model. 


Well  # 

Row,  Col 

Measured 

Head  (ft) 

Simulated 

Head  (ft) 

Measured  - 
Simulated 

Head  (ft) 

MP2MW1 

10,25 

679.69 

676.68 

3.01 

MP2MW3 

10,22 

675.51 

675.46 

0.05 

11,23 

675.34 

675.70 

-0.36 

MP2MW6 

1 1 ,20 

674.86 

674.52 

0.34 

CG3MW1 

6,24 

677.38 

676.69 

0.69 

CG3MW2 

7,22 

676.29 

675.78 

0.51 

CG3MW3 

9,23 

676.50 

675.98 

0.52 
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Well# 


CG3MW4_5* 


CG3MW7 


TF4MW1 


TF4MW2 


TF4MW3_4* 


SF5MW1 


SF5MW2 


SF5MW3_4* 


SF5MW6 


LF6MW1 


LF6MW2 


LF6MW3 


LF6MW4 


LF6MW5 


LF6MW6 


LF6MW8 


HN8MW1 


HN8MW2 


HN8MW3_4* 


RT9MW1 


RT9MW2 


RT9MW3 


RT9MW4_5* 


RT9MW6 


S1MW2 


S1MW3 


Row,  Col 


8,23 


8,20 


.9 


8,10 


9,10 


12,6 


14,5 


13,6 


13,5 


14,9 


14,8 


14,10 


16,7 


15,7 


13,10 


15,9 


5,22 


6,19 


7,21 


6,16 


7,14 


9,15 


8,14 


13,26 


13,27 


Measured 
Head  (ft) 


676.41 


675.64 


667.23 


659.61 


658.21 


674.15 


675.32 


676.46 


674.26 


672.1 


672.68 


671.93 


672.75 


673.07 


671.17 


673.12 


676.96 


675.31 


676.01 


673.06 


668.21 


670.72 


667.47 


670.58 


677.39 


677.15 


Simulated 

Head  (ft) 

Measured  • 
Simulated 

Head  (ft) 

676.08 

0.33 

674.80 

0.84 

668.83 

-1.60 

664.22 

-4.61 

660.06 

-1.85 

671.34 

2.81 

672.67 

2.65 

671.64 

4.82 

672.40 

1.86 

670.40 

1.70 

670.72 

1.96 

670.21 

1.72 

671.91 

0.84 

671.56 

1.51 

669.67 

1.50 

670.78 

2.34 

675.93 

1.03 

674.50 

0.81 

675.35 

0.66 

672.78 

0.28 

670.81 

-2.60 

671.26 

-0.54 

670.32 

-2,85 

672.33 

-1.75 

676.63 

0.76 

676.98 

0.17 
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Weil  « 

Row,  Col 

Measured 

Head  (ft) 

Simulated 

Head  (ft) 

Measured  • 
Simulated 

Head  (ft) 

S1MW11 

15,24 

675.72 

675.67 

0.05 

S1MW12 

16,25 

674.55 

675.94 

-1.39 

S1MW13 

14,24 

675.21 

675.76 

-0.55 

S1MW14 

14,25 

673.92 

676.14 

-2.22 

MP2MW2 

12,24 

675.57 

675.96 

-0.39 

*  Indicates  that  more  than  one  well  is  present  in  each  node  and  an  average  value  for 
hydraulic  head  was  used. 

Sum  of  Squared  Residuals  =  128.78/38  =  3.3891 
Root  Mean  Squared  Error  =  1.84 


It  should  also  be  noted  that  the  model  was  calibrated  with  respect  to  the  September  1993 
water  level  data  and  should  only  be  considered  calibrated  with  respect  to  this  data.  More 
information  on  the  water  level  fluctuation  through  time  would  be  needed  to  perform  a 
transient  calibration.  The  model  was  not  calibrated  with  respect  to  concentration  data,  but 
only  with  respect  to  hydraulic  head. 

Model  Assumptions  and  Limitations 

*  The  model  domain  consisted  only  of  the  shallow  unconfined  aquifer  (i.e.  one  layer). 

*  The  initial  head  data  input  to  the  transport  model  are  results  of  measurements  taken  in 
September  1993. 

*  Initial  concentrations  of  compounds  are  results  of  the  Round  IV  sampling  event  which 
was  conducted  from  July  to  September  1993. 

*  Hydraulic  conductivity  values  are  the  result  of  slug  tests  performed  in  November,  1987 
and  September  1993. 

*  The  model  was  calibrated  with  respect  to  hydraulic  head  using  September  1993  water 
level  data  and  should  only  be  considered  calibrated  with  respect  to  September  1993 
water  level  data. 

*  The  flow  model  was  assumed  to  be  at  steady-state  with  respect  to  hydraulic  head. 
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Site  2 


PCE  was  detected  in  two  wells  at  Site  2.  Wells  MP2MW7  and  MP2MW1  contained 
concentrations  of  PCE  of  6.3  and  7.2  ug/l,  respectively.  These  concentrations  were  input 
to  the  model  at  their  respective  nodes  and  the  concentrations  were  monitored  at 
Production  Well  #2  and  at  the  sinkhole  with  respect  to  time. 
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Appendix  R:  Site  3  Risk  Assessment 


iiiiiiiiiiiiiiiilllllilllliiiiri 


§i§§§§§§§§iiiiiiiiiiiii|||ili||| 

SSSSS5SSSdfNicNirsi(DN<bfMCD^d(NirM04<DW<Dr^CDr^OfMCNj 


^  0)  ^ 
<1-0 
T-  3  Cl 

O!  O  Ol 

O  D.O 

S 

ill  c  0) 

=  t 

E  < 

x  o 

n  7 


5  S  5 

f  f  ^ 

(0  (0  CO 
(0  CO  CO 

"  'o"  ^  'o" 

c  c  c 
0)  0)  a> 
OQ  OQ  CQ 


0)  0)  (U 

c  c  c 
0)  0)  0) 
o  o  o 

2  2  5 

.lit 

'  CO  CO  (0 

^  '.w*'  'w'  'w' 

o  o  o 

c  c  c 
a>  0)  0) 
CD  QQ  CO 


a> 

c 

0) 

£ 

0)  0  c 

C  C  CO 
0  0 

fe.  ^  s 

O.  Q.c^ 


O  O  O 
^  N  IsJ 
C  C  C 
0  0  0 
m  CD  CD 


0  0  0 
c  c  c 
0  0  0 

iff 

§  b  b 

<i.i.i 


o  o  o 

g  g  N 

0  0  0 
CQ  CD  CD 


r^oof^c^OO  OOOOOOO 

88S3SS3S  838888885g5588gg5g55S88« 


§888i888i88888|s|Sp|||S?S?||§|| 


SITE3R.WB1/ 940126 


CO 

£ 

CO 


o 

(/) 

JD 

3 

CO 

"O 

c 

CO 


O 

CO 


(/) 

£ 

,0 

2  i 

0)  CO 

^  o  c 
■o  E  0) 

'  o  -S- 


0) 


< 
*“0 


q: 

-S2 

CO 


3 

lo 

i2  S 

C  OJ 

=  f 

E  < 

CO 


CO 

c 

o 

(/) 

CO 

0) 

a: 


3 

CO 

O 


•o 

Q) 

N 


D 

CO 

CO 


!X 

H 

CL 

Ui  LU 

Q  CD 

*"1 

01 

< 

CO 


QQQQQooQ 0000000000 
000000000000000000 
920000000000000000 
ppoooooooooooooooo 
p  p  o  o’  d  o  o  d  d  d  d  d  d  d  d  d  d  d 


00000000000000 

00000000000000 

00000000000000 

00000000000000 

CO  d  d  d  d  d  d  d  d  d  iri  d  d  d 

r«^rv.r^t^h.r^h.r^r^incDr^r^r^ 


opopppoopooooooooooooooooooooooo 
oooooooooppooooooooooooooooooooo 
00000000000000000000000000000000 
ppppppppp 00000000000000000000000 
csi  <D  f\i  CD  cvi  (d  (>i  d  csi  <Ni  csi  (D  oj  CD  csi  CD  (N  d  fsi  fsi  cnI  (d  cNi  (o  c>j  CD  rsi  d  CN  rsi  fsi  CD 


99999000000000000000000000000000 
ooooopooopppopoooooooooooooooooo 
poopppoooooooooooooooooooooooooo 
popppppppooooooooooooooooooooooo 
d  d  Tt  d  Tf  d  (d  d  d  d  o W  d  d  00  d  o  d  d  d  ^  d  cd  d  d  d  ''t 


<D  0) 

c  c 

(U  (D 

(C  CO 


<D  0  0) 

C  C  C 

0)  0)  0) 

f  f  ^ 

(0  (0  (Q 

JL.  l_  k_ 

§  §  B 

^ 


a>  0)  0) 

c  c  c 

0)  0)  V 

(C  (0  (Q 

U.  L-  L- 

000 

3  3  3 


0)  0)  0) 

C  C  C 
<D  4)  O  _  , 

^  ^  ^  §J 

Q.  Q.  Q.  QJ 


O  O 
N  N 
C  C 

<U  4) 

CD  CD 


000 

N  N  N 
C  C  C 
4)  4)  4) 

CQ  03  CD 


000 

N  N  N 
C  C  C 

4;  4)  4) 

CD  CD  CD 


4;  4)  4) 
C  C  C 
4)  4)  4) 
(/)  (0  (0 
>S  >N  >s 


C 

CO 

^  a 

c  o 

4;  N 

w  c 

k.  X 

x: 


x:  x:  ^  ^  X  ^  ^  X  X  _  _ 

OOOOOOOOOOQQQQQQCI^CIiQ 


T3 

o” 

^  rsj 

4)  4) 

-S 

Q  c 


*D  *3  -D 
o"  o'  o" 

■  I  I 

CO  CO  CO 
CN  CnT  CvT 


000 
c  c  c 

4}  O  4) 
T3  -3  -D 
C  C  C 


0000000  0000000  0000000 

CDCDCOCOCOCDCO  ,  ,  .COCOCOCOCOCOCO  .  .  .COCOCOCOCOCOCO  CD 

^CDCDCDCDCDCDCD^z^=dcDCQCaCDCDCaCD^d=!cDCDCDCDCDCDCD=^d=iCD 

OD3DDD3D00ODDZ)3DDDOOO3DDDDDDO0O3 

CO(O(/)C/)(DCOCOC/)CO(DC0COC/}(O  COC/)CO(/)COCOCOCOC/)COC/3C/)tOC/)C/)(/)C/)C/> 


(NCDCNCDCNICOCNO 

00'r-0'*-0't-T- 
I  I  I  ■  I  I  I  I 

O^OTfOTfOOO 

OOt-Ot-O'^-O 

(OCDCOCOCOCOCOCO 

COCOCOCOCDCDCDCO 


r^CNJCNCDCNJCOCNCDCgO 

OOOOt-O’i-O’t-t- 


(NCsirMCDCMCOCNJCDCNOCNfMOaCD 

0000'»-Ot“0'»-t-0000 


oaoTj-oxfO"^ 

0000'^0'«-0 

COCDCOCOCDCOCOCO 

COCDCOCOCOCDCOCO 


CO 


T-  <N  (N  CO  CO  CO 


CD  CD  CD  m 
CD  CD  CD  CO 

till 

CO  CO  CO  CO 

CD  CD  CD  CD 

O  O  O  U 

■  ill 

0000 

a  CL  CL  CL 


CD  CD  CD  CD 
CD  CD  CD  CD 

I  I  I  I 

CO  CO  CO  CO 

CD  CD  CD  CD 

0000 

6  6  6  6 

a  Q-  Q.  CL 


CNJ  CO  "r- 

m  m  m  m 

CD  CD  CD  CD 

I  I  I  I 

CO  CO  CO  CO 

0  CD  CD  CD 

0000 

6  6  6  6 

CL  Q.  CL  a 


oooooo'^O'^OTfocoood'^ 
t-OOOOOT-OT-O-t-OOOOO 
CDCD00(D00CDW0CD00CDW0 
CDCDCDCDCDWppWCDCDCDCDCDWCD 

t^CNJCNCOCOCOt^CNCO’ . 


CN  CN  CO  CO  CO 


CM  CO 


CD  CD  CD  CD 
0  CD  CD  0 

I  I  I  I 

CO  CO  0  CO 
0  0  0  0 
0000 

till 

0000 

CL  CL  Q.  CL 


0  CD  0  0  CD 

0  0  0  0  0 

I  I  I  I  I 
0  0  0  0  0 

0  0  0  0  0 

00000 

6  6  6  6  6 

CL  CL  CL  CL  Q_ 


0  0  0  0  0 

CD  CD  CD  CD  CD 

0  0  0  0  CO 

0  0  0  0  0 

00000 

I  I  I  I  I 

00000 

CL  CL  CL  CL  Q. 


0  0  0  0  0  0 
CD  CD  CD  CD  CD  CD 
0  0  0  0  0  0 
0  0  0  0  0  0 
000000 

I  I  I  t  I  I 

000000 

CL  CL  CL  a  CL  CL 


0'r-«f-fNJCNI000T-rsJ0 


T-(Nrsj000T-rsj0 


Cvl  CM  0  0  0 


CM  0 


000000000000000000000300000000000 
000000000000  00000000000000000000 


CDOOCDO(DOCDCDCDCDCDCDCDCDCD(DCDCDCDCD0CDCD0CDCD0CDCD(DCD 

00000000000000000000000000000000 

00000000000000000000000000000000 

Q.  Q.  Q-  Q.  Q.  CLQ.Q.Q.|]LQ.Q.Q.Q.Q.a.CLQ.Q.CLQ.Q.Q.£LQ.aa.Q.Q.Q.Q.Q. 


R-2 


SITE3R.WB1/ 940126 


SltE3R.WB1/ 940126 


TABLE  R-1  EXPOSURE  ASSESSMENT  PARAMETERS  -  DERMAL  CONTACT  WITH  SOILS 
Site  3,  Former  County  Garage,  Alpena  CRTC,  Alpena,  Ml 
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TABLE  R-4  EXPOSURE  ASSESSMENT  PARAMETERS  -  SOIL  INGESTION 
Site  3.  Former  County  Garage,  Alpena  CRTC,  Alpena,  Ml 
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site  3,  Fortner  County  Garage,  Alpena  CRTC,  Alpena,  Ml 


TABLE  B-11  SOIL  INHALATION  MODEL  CALCULATION 

SK«  3,  Eorm«r  County  Oirage,  Alpana  CRTC,  Alpona,  Ml 
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Groundwater  Solute  Transport  Model  Data  -  Site  3 
Theoretical  Background 

A  two-dimensional  Method  of  Characteristics  (MOC)  solute  transport  model  (Konikow  and 
Bredehoeft,  1989)  was  used  for  preliminary  examination  of  contaminant  migration  within 
the  shallow  aquifer  beneath  the  Alpena  CRTC.  The  model  is  designed  to  calculate 
transient  changes  in  solute  concentrations  within  groundwater  by  simultaneously  solving 
partial  differential  equations  describing  groundwater  flow  and  transport  and  computes  the 
change  in  a  chemicals  concentration  over  time.  Changes  in  chemical  concentrations  over 
time  are  caused  by  the  processes  of  convective  transport,  hydrodynamic  dispersion,  and 
mixing  from  fluid  sources.  This  model  couples  the  groundwater  flow  equation  with  solute 
transport  equations. 

The  flow  equation  can  be  approximated  by  an  implicit  finite-difference  equation.  The 
model  area  is  discretized  into  a  rectangular  grid  with  each  square  being  a  node.  The  finite 
difference  equation  is  solved  numerically  for  each  node  in  the  grid  using  an  iterative 
alternating-direction  implicit  (ADI)  procedure. 

After  the  hydraulic  head  distribution  is  calculated,  the  velocity  of  groundwater  flow  can  be 
computed  at  each  node.  The  expression  for  average  velocity  of  groundwater  can  be 
derived  from  Darcy's  law.  The  groundwater  velocity  at  each  node  is  calculated  utilizing  an 
explicit  finite-difference  approximation  of  Darcy's  law.  The  computer  program  uses  an 
alternating-direction  implicit  procedure  to  solve  a  finite-difference  approximation  to  the 
groundwater  flow  equation,  and  it  uses  the  method  of  characteristics  (MOC)  to  solve  the 
solute  transport  equation.  MOC  uses  a  particle  tracking  procedure  to  represent  convective 
transport  and  a  two-step,  explicit  procedure  to  solve  a  finite-difference  equation  that 
describes  the  effects  of  hydrodynamic  dispersion,  fluid  sources  and  sinks,  and  divergence 
of  velocity. 

A  number  of  assumptions  are  inherent  in  the  solute  transport  model: 

1 .  Darcy's  law  is  valid  and  hydraulic  head  gradients  are  the  only  significant  driving 
mechanism  for  fluid  flow. 

2.  The  porosity  and  hydraulic  conductivity  of  the  aquifer  are  constant  with  time,  and 
porosity  is  uniform  in  space. 

3.  Gradients  of  fluid  density,  viscosity,  and  temperature  do  not  affect  the  velocity 
distribution. 

4.  No  chemical  reactions  occur  that  affect  the  concentration  of  the  solute,  the  fluid 
properties,  or  the  aquifer  properties. 

5.  Ionic  and  molecular  diffusion  are  negligible  contributors  to  the  total  dispersive  flux. 
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6.  Vertical  variations  in  head  and  concentrations  are  negligible. 

7.  The  aquifer  is  homogeneous  and  isotropic  with  respect  to  the  coefficients  of 
longitudinal  and  transverse  dispersivity. 


Transport  Model  Input 

A  model  grid  of  32  columns  by  19  rows  with  a  250  foot  lateral  spacing  was  used. 
Specified  head  cells  were  used  at  nodes  corresponding  to  the  South  Branch  of  the  Thunder 
Bay  River,  at  nodes  along  the  eastern  boundary  of  the  model  grid  area,  and  also  at  the 
sinkhole  in  the  northeastern  portion  of  the  model.  Groundwater  elevations  measured 
during  September,  1993  were  used  as  initial  input  into  the  transport  model. 

Hydraulic  conductivity  values  were  calculated  from  slug  tests  performed  at  Alpena  CRTC 
(Engineering  Science,  1989;  Earth  Technology,  1994).  Values  of  hydraulic  conductivity 
range  from  1 2  feet/day  at  Site  4  to  278  feet/day  at  Site  3. 

Aquifer  thickness  values  were  obtained  from  drilling  records  of  monitoring  wells  and  soil 
borings  obtained  from  the  SI  and  Rl  field  activities.  Values  listed  are  from  logs  in  which 
the  thickness  of  the  shallow  aquifer  was  clearly  discernible,  and  ranged  from  20  feet  at 
Site  5  to  65  feet  at  Site  8.  Transmissivity  values  were  calculated  by  multiplying  the 
calculated  hydraulic  conductivity  values  by  the  aquifer  thickness.  Transmissivity  ranges 
from  420  ft^/day  at  TF4-MW3  to  1  5,290  ft^/day  at  CG3-MW5. 

Monitoring  of  the  discharge  of  springs  into  the  sinkhole  was  performed  during  the  SI 
(Engineering  Science,  1990)  and  an  estimate  of  approximately  18,000  gallons  of  water  per 
day  discharging  into  the  sinkhole  was  calculated.  In  order  to  obtain  a  numerical  estimate 
of  discharge  into  the  sinkhole  for  the  model,  MODFLOW  (McDonald  and  Harbaugh,  1988), 
a  3-  dimensional  finite-difference  groundwater  flow  model  was  used.  MODFLOW  was 
used  because  of  its  ability  to  simulate  the  effect  of  head-dependent  groundwater  flow  into 
a  groundwater  sink  (i.e.  the  sinkhole).  This  package  was  not  available  in  MOC.  The  same 
model  parameters  and  boundary  conditions  were  used  within  MODFLOW  as  in  MOC. 

Based  upon  hydraulic  head  data  collected  in  September  1993,  discharge  from  the  shallow 
aquifer  into  the  sinkhole  is  approximately  30,000  gallons  per  day. 

The  dispersivity  of  an  aquifer  in  two  dimensions  is  described  by  the  longitudinal  dispersion, 
the  transverse  dispersion  and  the  ratio  of  the  two  (Fetter,  1993).  As  a  contaminant  plume 
moves  further  from  its  initial  location  within  the  aquifer  by  advection  with  the  groundwater 
flow,  the  plume  spreads.  The  spreading  in  the  direction  of  groundwater  flow  is  the 
longitudinal  dispersion,  the  spreading  in  the  direction  perpendicular  to  the  groundwater 
flow  is  known  as  the  transverse  dispersion  (Fetter,  1993).  The  values  of  the  dispersion 
coefficients  are  typically  derived  via  bench  scale  tests,  aquifer  tests,  or  calibration  of 
contaminant  transport  models.  Since  no  data  presently  exists  describing  dispersivity 
within  the  shallow  aquifer  beneath  the  Alpena  CRTC  and  insufficient  data  exists  to  allow 
for  derivation  of  dispersivity  via  model  calibration,  moderate  values  of  100  feet  for 
longitudinal  dispersivity  and  30  feet  for  transverse  dispersivity  were  chosen  (Gillham  and 
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Cherry,  1 982).  A  more  complete  description  of  the  model  is  given  in  the  report. 
Preliminary  Groundwater  Modeling  Effort,  Earth  Technology,  August  1993. 


Model  Calibration: 

The  groundwater  flow  model  was  calibrated  with  respect  to  the  September  1993 
groundwater  elevation  measurements.  Calibration  of  the  groundwater  flow  model  was 
accomplished  by  defining  a  set  of  parameters,  boundary  conditions,  and  stresses  that 
produce  simulated  heads  and  fluxes  that  match  field-measured  values  within  a 
preestablished  range  of  error  (Anderson  and  Woessner,  1992).  In  order  to  match  field 
measured  values  for  hydraulic  head  as  determined  during  September  1993,  a  few 
modifications  were  made  to  the  preexisting  groundwater  flow  model.  These  changes 
included  updating  the  initial  head  array,  modeling  the  sinkhole  as  a  constant  head  cell  to 
account  for  the  large  gradient  changes  in  the  vicinity  of  the  sinkhole  and  including  recharge 
to  the  model  at  a  rate  of  9  inches  per  year  over  the  whole  model  area.  By  adjusting  these 
parameters,  an  acceptable  level  of  calibration  was  achieved.  An  acceptable  level  of 
calibration  was  defined  as  a  root  mean  squared  error  (RMS)  of  less  than  2  feet.  The  RMS, 
or  the  standard  deviation  is  the  average  squared  difference  in  measured  and  simulated 
heads  and  is  given  by  the  equation: 


RMS=[l/nY, 

i-l 


n  =  number  of  wells 
hm  =  measured  head 
h,  =  model  simulated  head 


Certain  portions  of  the  model  may  have  values  above  the  goal  of  2  feet  while  others  fall 
much  below  this  value.  The  RMS  represents  the  average  error  present  in  the  model. 

The  following  provides  a  summary  of  the  final  calibrated  heads  for  the  flow  model. 


Well  # 

Row,  Col 

Measured 

Head  (ft) 

Simulated 

Head  (ft) 

Measured  - 
Simulated 

Head  (ft) 

MP2MW1 

10,25 

679.69 

676.68 

3.01 

MP2MW3 

10,22 

675.51 

675.46 

0.05 

MP2MW4_5* 

11,23 

675.34 

675.70 

-0.36 

MP2MW6 

1 1 ,20 

674.86 

674.52 

0.34 

CG3MW1 

6,24 

677.38 

676.69 

0.69 

CG3MW2 

7,22 

676.29 

675.78 

0.51 
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Wen# 

Row,  Col 

Measured 

Head  (ft) 

Simulated 

Head  (ft) 

Measured  • 
Simulated 

Head  (ft) 

CG3MW3 

9,23 

676.50 

675.98 

0.52 

CG3MW4_5* 

8,23 

676.41 

676.08 

0.33 

CG3MW7 

8,20 

675.64 

674.80 

0.84 

TF4MW1 

7,9 

667.23 

668.83 

-1.60 

TF4MW2 

8,10 

659.61 

664.22 

-4.61 

TF4MW3_4* 

9,10 

658.21 

660.06 

-1.85 

SF5MW1 

12,6 

674.15 

671.34 

2.81 

SF5MW2 

14,5 

675.32 

672.67 

2.65 

SF5MW3_4* 

13,6 

676.46 

671.64 

4.82 

SF5MW6 

13,5 

674.26 

672.40 

1.86 

LF6MW1 

14,9 

672.1 

670.40 

1.70 

LF6MW2 

14,8 

672.68 

670.72 

1.96 

LF6MW3 

14,10 

671.93 

670.21 

1.72 

LF6MW4 

16,7 

672.75 

671.91 

0.84 

LF6MW5 

15,7 

673.07 

671.56 

1.51 

LF6MW6 

13,10 

671.17 

669.67 

1.50 

LF6MW8 

15,9 

673.12 

670.78 

2.34 

HN8MW1 

5,22 

676.96 

675.93 

1.03 

HN8MW2 

6,19 

675.31 

674.50 

0.81 

HN8MW3_4* 

7,21 

676.01 

675.35 

0.66 

RT9MW1 

6,16 

^  673.06 

672.78 

0.28 

RT9MW2 

7,14 

668.21 

670.81 

-2.60 

RT9MW3 

9,15 

670.72 

671.26 

-0.54 

RT9MW4_5* 

8,14 

667.47 

670.32 

-2.85 

RT9MW6 

8,16  - 

670.58 

672.33 

-1.75 

S1MW2 

13,26 

677.39 

676.63 

0.76 
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*  Indicates  that  more  than  one  well  is  present  in  each  node  and  an  average  value  for 
hydraulic  head  was  used. 

Sum  of  Squared  Residuals  =  128.78/38  =  3.3891 
Root  Mean  Squared  Error  =  1 .84 


It  should  also  be  noted  that  the  model  was  calibrated  with  respect  to  the  September  1993 
water  level  data  and  should  only  be  considered  calibrated  with  respect  to  this  data.  More 
information  on  the  water  level  fluctuation  through  time  would  be  needed  to  perforrn  a 
transient  calibration.  The  model  was  not  calibrated  with  respect  to  concentration  data,  bu 
only  with  respect  to  hydraulic  head. 

Model  Assumptions  and  Limitations 

*  The  model  domain  consisted  only  of  the  shallow  unconfined  aquifer  (i.e.  one  layer). 

*  The  initial  head  data  input  to  the  transport  model  are  results  of  measurements  taken  in 
September  1 993. 

*  Initial  concentrations  of  compounds  are  results  of  the  Round  IV  sampling  event  which 
was  conducted  from  July  to  September  1 993. 

*  Hydraulic  conductivity  values  are  the  result  of  slug  tests  performed  in  November.  1 987 
and  September  1993. 

*  The  model  was  calibrated  with  respect  to  hydraulic  head  using  September  1 993  water 
level  data  and  should  only  tie  considered  calibrated  with  respect  to  September  1993 
water  level  data. 

*  The  flow  model  was  assumed  to  be  at  steady-state  with  respect  to  hydraulic  head. 
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Site  3 


No  organic  compounds  or  inorganics  exceeding  the  Act  307  Type  A  or  B  cleanup  criteria 
were  present  in  the  shallow  aquifer  at  this  site. 
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Appendix  S:  Site  4  Risk  Assessment 


Table  S-1A 

a  Utilized  to  Calculate  Reasonable  Maximum  Exposure  Concentrations  for  Groundwater  -  Site  4 

MIANG,  Alpena  CRTC,  Alpena,  Ml 
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Table  S-2A 

Data  Utilized  to  Calculate  Reasonable  Maximum  Exposure  Concentrations  for  Sediment  -  Site  4 

MIANG,  Alpena  CRTC,  Alpena,  Ml 
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TABLE  S-1  EXPOSURE  ASSESSMENT  PARAMETERS  -  DERMAL  CONTACT  WITH  SEDIMENTS 
SHb  4,  Third  Fire  Training  Area,  Alpena  CRTC,  Alpena,  Ml 
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TABLE  S-4  DAILY  INTAKE  -  ADULT 

DERMAL  CONTACT  WITH  SEDIMENTS 
FUTURE  LAND  USE  SCENARIO 
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TABLE  S-21  EXPOSURE  ASSESSMENT  PARAMETERS  -  DERMAL  CONTACT  WITH  THUNDERBAY  RIVER  SURFACE  WATER 
Site  4,  Third  Fire  Training  Area,  Alpena  CRTC,  Alpena,  Ml 
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TABLE  S-24  DAILY  INTAKE  -  ADULT 

SURFACE  WATER  DERMAL  CONTACT 
FUTURE  LAND  USE  SCENARIO 
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TABLE  S-27  SUBCHRONK:  HAZARDOUS  INDEX  ESTWATE  -  DERMAL  CONTACT  WITH  SEDIMENTS 
FUTURE  LAND  USE  SCENARIO  -  ADULTS/CHILDREN 


TABLE  S-28  SUBCHRONIC  HAZARDOUS  INDEX  ESTIMATE  -  SEDIMENT  INGESTION 
FUTURE  LAND  USE  SCENARIO  •  ADULTS/CHILDREN 
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TABLE  S-31  CHRONIC  HAZARDOUS  INDEX  ESTIMATE  -  FISH  INGESTION 
FUTURE  LAND  USE  SCENARIO  -  ADULTS  AND  CHILDREN 
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TABLE  5-33  CHRONIC  HAZARDOUS  INDEX  ESTIMATE  -  SURFACE  WATER  INGESTION 
FUTURE  LAND  USE  SCENARIO  -  ADULTS  AND  CHILDREN 
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TABLE  S-35  CHRONIC  HAZARDOUS  INDEX  ESTIMATE  •  SURFACE  WATER  DERMAL  CONTACT 
FUTURE  LAND  USE  SCENARIO  -  ADULTS  AND  CHILDREN 
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Groundwater  Solute  Transport  Model  Data  -  Site  4 


Theoretical  Background 

A  two-dimensional  Method  of  Characteristics  (MOO  solute  transport  model  (Konikow  and 
Bredehoeft,  1989)  was  used  for  preliminary  examination  of  contaminant  migrahon  with 
The  shallow  aquifer  beneath  the  Alpena  CRTC.  The  model  is  designed  to  calculate 
transient  changes  in  solute  concentrations  within  groundwater  by  simultaneously  solving 
oartfal  differen?ial  equations  describing  groundwater  flow  and  transport  and  computes  the 
Ke  in  a  c^mic^s  concentration  Ivor  time.  Changes  in  chemical  concentrations  over 
time  are  caused  by  the  processes  of  convective  transport,  hydrodynamic  dispersion,  and 
mixing  from  fluid  sources.  This  model  couples  the  groundwater  flow  equation  with  solute 

transport  equations. 

The  flow  equation  can  be  approximated  by  an  implicit  finite-difference  equation.  Jhe 
mr^al  a^aTs  dlLratizad  into  a  rectangular  grid  vrrith  each 

difference  equation  is  solved  numerically  for  each  node  in  the  grid  using  an  iterative 
alternating-direction  implicit  (ADI)  procedure. 

After  the  hydraulic  head  distribution  is  calculated,  the  velocity  of  groundwater  flow  can  be 
romnuted  at  each  node.  The  expression  for  average  velocity  of  groundwater  can  be 
derived  from  Darcy's  law.  The  groundwater  velocity  at  each  node  is  calculated  utilizing  an 
erpncit  finTe-difference  approximation  of  Darcy's  law.  The  computer 
alternating-direction  implicit  procedure  to  solve  a  finite-difference  °  . 

oroundwater  flow  equation,  and  it  uses  the  method  of  characteristics  (MOC)  to  solve  the 
S  tta  “sp«"quron.  MOC  uses  a  particle  tracking  procedure  to  represent  convective 
transport  and  a  two-step,  explicit  procedure  to  solve  a  finite-difference  equation  that 
describes  the  effects  of  hydrodynamic  dispersion,  fluid  sources  and  sinks,  and  diverge 

of  velocity. 

A  number  of  assumptions  are  inherent  in  the  solute  transport  model: 

1 .  Darcy's  law  is  valid  and  hydraulic  head  gradients  are  the  only  significant  driving 
mechanism  for  fluid  flow. 

2.  The  porosity  and  hydraulic  conductivity  of  the  aquifer  are  constant  with  time,  and 
porosity  is  uniform  in  space. 

3.  Gradients  of  fluid  density,  viscosity,  and  temperature  do  not  affect  the  velocity 
distribution. 


4.  No  chemical  reactions  occur  that  affect  the  concentration  of  the  solute,  the  fluid 
properties,  or  the  aquifer  properties. 

5.  Ionic  and  molecular  diffusion  are  negligible  contributors  to  the  total  dispersive  flux. 
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6.  Vertical  variations  in  head  and  concentrations  are  negligible. 

7.  The  aquifer  is  homogeneous  and  isotropic  with  respect  to  the  coefficients  of 
longitudinal  and  transverse  dispersivity. 


Transport  Model  Input 

A  model  grid  of  32  columns  by  1 9  rows  with  a  250  foot  lateral  spacing  was  used. 
Specified  head  cells  were  used  at  nodes  corresponding  to  the  South  Branch  of  the  Thunder 
Bay  River,  at  nodes  along  the  eastern  boundary  of  the  model  grid  area,  and  also  at  the 
sinkhole  in  the  northeastern  portion  of  the  model.  Groundwater  elevations  measured 
during  September,  1 993  were  used  as  initial  input  into  the  transport  model. 

Hydraulic  conductivity  values  were  calculated  from  slug  tests  performed  at  Alpena  CRTC 
(Engineering  Science,  1989;  Earth  Technology,  1994).  Values  of  hydraulic  conductivity 
range  from  1 2  feet/day  at  Site  4  to  278  feet/day  at  Site  3. 

Aquifer  thickness  values  were  obtained  from  drilling  records  of  monitoring  wells  and  soil 
borings  obtained  from  the  SI  and  Rl  field  activities.  Values  listed  are  from  logs  in  which 
the  thickness  of  the  shallow  aquifer  was  clearly  discernible,  and  ranged  from  20  feet  at 
Site  5  to  65  feet  at  Site  8.  Transmissivity  values  were  calculated  by  multiplying  the 
calculated  hydraulic  conductivity  values  by  the  aquifer  thickness.  Transmissivity  ranges 
from  420  ft^/day  at  TF4-MW3  to  15,290  ftVday  at  CG3-MW5. 

Monitoring  of  the  discharge  of  springs  into  the  sinkhole  was  performed  during  the  SI 
(Engineering  Science,  1990)  and  an  estimate  of  approximately  18,000  gallons  of  water  per 
day  discharging  into  the  sinkhole  was  calculated.  In  order  to  obtain  a  numerical  estimate 
of  discharge  into  the  sinkhole  for  the  model,  MODFLOW  (McDonald  and  Harbaugh,  1988), 
a  3-  dimensional  finite-difference  groundwater  flow  model  was  used.  MODFLOW  was 
used  because  of  its  ability  to  simulate  the  effect  of  head-dependent  groundwater  flow  into 
a  groundwater  sink  (i.e.  the  sinkhole).  This  package  was  not  available  in  MOC.  The  same 
model  parameters  and  boundary  conditions  were  used  within  MODFLOW  as  in  MOC. 

Based  upon  hydraulic  head  data  collected  in  September  1993,  discharge  from  the  shallow 
aquifer  into  the  sinkhole  is  approximately  30,000  gallons  per  day. 

The  dispersivity  of  an  aquifer  in  two  dimensions  is  described  by  the  longitudinal  dispersion, 
the  transverse  dispersion  and  the  ratio  of  the  two  (Fetter,  1993).  As  a  contaminant  plume 
moves  further  from  its  initial  location  within  the  aquifer  by  advection  with  the  groundwater 
flow,  the  plume  spreads.  The  spreading  in  the  direction  of  groundwater  flow  is  the 
longitudinal  dispersion,  the  spreading  in  the  direction  perpendicular  to  the  groundwater 
flow  is  known  as  the  transverse  dispersion  (Fetter,  1993).  The  values  of  the  dispersion 
coefficients  are  typically  derived  via  bench  scale  tests,  aquifer  tests,  or  calibration  of 
contaminant  transport  models.  Since  no  data  presently  exists  describing  dispersivity 
within  the  shallow  aquifer  beneath  the  Alpena  CRTC  and  insufficient  data  exists  to  allow 
for  derivation  of  dispersivity  via  model  calibration,  moderate  values  of  100  feet  for 
longitudinal  dispersivity  and  30  feet  for  transverse  dispersivity  were  chosen  (Gillham  and 
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Cherry,  1982).  A  more  complete  description  of  the  model  is  given  in  the  report. 
Preliminary  Groundwater  Modeling  Effort,  Earth  Technology,  August  1993. 


Model  Calibration: 

The  groundwater  flow  model  was  calibrated  with  respect  to  the  September  1993 
groundwater  elevation  measurements.  Calibration  of  the  groundwater  flow  model  was 
accomplished  by  defining  a  set  of  parameters,  boundary  conditions,  and  stresses  that 
produce  simulated  heads  and  fluxes  that  match  field-measured  values  within  a 
preestablished  range  of  error  {Anderson  and  Woessner,  1992).  In  order  to  match  field 
measured  values  for  hydraulic  head  as  determined  during  September  1 993,  a  few 
modifications  were  made  to  the  preexisting  groundwater  flow  model.  These  changes 
included  updating  the  initial  head  array,  modeling  the  sinkhole  as  a  constant  head  cell  to 
account  for  the  large  gradient  changes  in  the  vicinity  of  the  sinkhole  and  including  recharge 
to  the  model  at  a  rate  of  9  inches  per  year  over  the  whole  model  area.  By  adjusting  these 
parameters,  an  acceptable  level  of  calibration  was  achieved.  An  acceptable  level  of 
calibration  was  defined  as  a  root  mean  squared  error  (RMS)  of  less  than  2  feet.  The  RMS, 
or  the  standard  deviation  is  the  average  squared  difference  in  measured  and  simulated 
heads  and  is  given  by  the  equation: 


RMS=[l/n'^ 

1*1 


n  =  number  of  wells 
h„  =  measured  head 
h,  =  model  simulated  head 

Certain  portions  of  the  model  may  have  values  above  the  goal  of  2  feet  while  others  fall 
much  below  this  value.  The  RMS  represents  the  average  error  present  in  the  model. 

The  following  provides  a  summary  of  the  final  calibrated  heads  for  the  flow  model. 


Well  # 

Row,  Col 

Measured 

Head  (ft) 

Simulated 

Head  (ft) 

Measured  - 
Simulated 

Head  (ft) 

MP2MW1 

10,25 

679.69 

676.68 

3.01 

MP2MW3 

10,22 

675.51 

675.46 

0.05 

11,23 

675.34 

675.70 

-0.36 

MP2MW6 

11,20 

674.86 

674.52 

0.34 

CG3MW1 

6,24 

677.38 

676.69 

0.69 

CG3MW2 

7,22 

676.29 

675.78 

0.51 
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Well  # 


CG3MW3 


CG3MW4_5* 


CG3MW7 


TF4MW1 

TF4MW2 

TF4MW3_4« 

SF5MW1 

SF5MW2 

SF5MW3_4* 

SF5MW6 

LF6MW1 

LF6MW2 

LF6MW3 

LF6MW4 

LF6MW5 

LF6MW6 

LF6MW8 


HN8MW1 


HN8MW2 


HN8MW3_4* 


RT9MW1 


RT9MW2 


RT9MW3 


RT9MW4_5* 


RT9MW6 


S1MW2 


Row,  Col 


9,23 


8,23 


8,20 


,9 


8,10 


9,10 


12,6 


14,5 


13,6 


13,5 


14,9 


14,8 


14,10 


16,7 


15,7 


13,10 


15,9 


5,22 


6,19 


7,21 


6,16 


7,14 


9,15 


8,14 


6 


13,26 


Measured 

Head  (ft) 

Simulated 

Head  (ft) 

676.50 

675.98 

676.41 

676.08 

675.64 

674.80 

667.23 

668.83 

659.61 

664.22 

658.21 

660.06 

674.15 

671.34 

675.32 

672.67 

676.46 

671.64 

674.26 

672.40 

672.1 

670.40 

672.68 

670.72 

671.93 

670.21 

672.75 

671.91 

673.07 

671.56 

671.17 

669.67 

673.12 

670.78 

676.96 

675.93 

675.31 

674.50 

676.01 

675.35 

673.06 

672.78 

668.21 

670.81 

670.72 

671.26 

667.47 

670.32 

670.58 

672.33 

677.39 

676.63 

0.33 


0.84 


-1.60 


-4.61 


-1.85 


2.81 


2.65 


4.82 


1.86 


1.70 


.96 


1.72 


0.84 


1.51 


1.50 


1.03 


0.81 


0.66 


0.28 


-2.60 


-0.54 


-2.85 


-1.75 


Weil# 

Row,  Col 

Measured 

Head  (ft) 

Simulatjed 

Head  (ft) 

Measured  - 
Simulated 

Head  (ft) 

S1MW3 

13,27 

677.15 

676.98 

0.17 

S1MW11 

15,24 

675.72 

675.67 

0.05 

S1MW12 

16,25 

674.55 

675.94 

-1.39 

S1MW13 

14,24 

675.21 

675.76 

-0.55 

S1MW14 

14,25 

673.92 

676.14 

-2.22 

MP2MW2 

12,24 

675.57 

675.96 

-0.39 

*  Indicates  that  more  than  one  well  is  present  in  each  node  and  an  average  value  for 
hydraulic  head  was  used. 

Sum  of  Squared  Residuals  =  128.78/38  =  3.3891 
Root  Mean  Squared  Error  =  1 .84 


It  should  also  be  noted  that  the  model  was  calibrated  with  respect  to  the  September  1993 
water  level  data  and  should  only  be  considered  calibrated  with  respect  to  this  data.  More 
information  on  the  water  level  fluctuation  through  time  would  be  needed  to  perform  a 
transient  calibration.  The  model  was  not  calibrated  with  respect  to  concentration  data,  but 
only  with  respect  to  hydraulic  head. 

Model  Assumptions  and  Limitations 

*  The  model  domain  consisted  only  of  the  shallow  unconfined  aquifer  (i.e.  one  layer). 

*  The  initial  head  data  input  to  the  transport  model  are  results  of  measurements  taken  in 
September  1993. 

*  Initial  concentrations  of  compounds  are  results  of  the  Round  IV  sampling  event  which 
was  conducted  from  July  to  September  1993. 

*  Hydraulic  conductivity  values  are  the  result  of  slug  tests  performed  in  November,  1987 
and  September  1993. 

*  The  model  was  calibrated  with  respect  to  hydraulic  head  using  September  1993  water 
level  data  and  should  only  be  considered  calibrated  with  respect  to  September  1993 
water  level  data. 

*  The  flow  model  was  assumed  to  be  at  steady-state  with  respect  to  hydraulic  head. 
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Site  4 


No  organic  compounds  were  present  in  the  groundwater  at  Site  4  above  Act  307  Type  B 
cleanup  criteria.  However,  model  simulations  predict  that  contaminants  from  sites  2,  5,  6, 
8,  and  9  would  reach  the  sinkhole  after  a  period  of  several  years  with  maximum 
concentrations  entering  the  sinkhole  between  5  and  20  years  depending  on  which  site  the 
contaminants  originate  from.  PCE,  arsenic,  and  carbon  tetrachloride  originating  from  wells 
MP2MW7,  MP2MW1  and  PW3  at  Site  2  arrive  at  the  sinkhole  in  concentrations  ranging 
from  0.002  (arsenic)  to  0.28  ug/I  for  PCE.  The  chemicals  1 ,2  -DCA,  1 ,4  -DCB,  benzene 
and  styrene  originating  from  well  SF5MW1  at  Site  5  arrive  at  the  sinkhole  at 
concentrations  ranging  from  0.01  (styrene)  and  1.55  ug/I  for  benzene.  Well  LF6MW3  at 
Site  6  contained  carbon  tetrachloride  which  arrived  at  the  sinkhole  at  an  estimated 
concentration  of  0.08  ug/I.  Well  HN8MW3  at  Site  8  contained  PCE  which  arrived  at  the 
sinkhole  with  an  estimated  concentration  of  0.001  ug/I.  Well  RT9MW6  at  Site  9 
contained  PCE,  benzene,  1,4  -DCB,  lead,  and  2-methyInapthalene  which  arrived  at  the 
sinkhole  at  concentrations  ranging  from  0.001  (PCE)  to  0.01  ug/I  for  2-methylnapthlene. 
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Appendix  T:  Site  5  Risk  Assessment 
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EXPOSURE  ASSESSMENT  PARAMETERS  -  FISH  INGESTION 
e  6,  Second  Fire  Training  Area,  Alpena  CRTC,  Alpena,  Ml 
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Carcinogenic  Effect  Exposure  Assessment 

Site  5,  Second  Fire  Training  Area,  Alpena  CRTC,  Alpena,  Ml 
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TABLE  T-7  CHRONIC  HAZARDOUS  INDEX  ESTIMATE  -  FISH  INGESTION 
FUTURE  LAND  USE  SCENARIO  •  ADULTS  AND  CHILDREN 
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Groundwater  Solute  Transport  Model  Data  -  Site  5 
Theoretical  Background 

A  two-dimensional  Method  of  Characteristics  (MOC)  solute  transport  model  (Konikow  and 
Bredehoeft,  1989)  was  used  for  preliminary  examination  of  contaminant  migration  within 
the  shallow  aquifer  beneath  the  Alpena  CRTC.  The  model  is  designed  to  calculate 
transient  changes  in  solute  concentrations  within  groundwater  by  simultaneously  solving 
partial  differential  equations  describing  groundwater  flow  and  transport  and  computes  the 
change  in  a  chemicals  concentration  over  time.  Changes  in  chemical  concentrations  over 
time  are  caused  by  the  processes  of  convective  transport,  hydrodynamic  dispersion,  and 
mixing  from  fluid  sources.  This  model  couples  the  groundwater  flow  equation  with  solute 
transport  equations. 

The  flow  equation  can  be  approximated  by  an  implicit  finite-difference  equation.  The 
model  area  is  discretized  into  a  rectangular  grid  with  each  square  being  a  node.  The  finite 
difference  equation  is  solved  numerically  for  each  node  in  the  grid  using  an  iterative 
alternating-direction  implicit  (ADI)  procedure. 

After  the  hydraulic  head  distribution  is  calculated,  the  velocity  of  groundwater  flow  can  be 
computed  at  each  node.  The  expression  for  average  velocity  of  groundwater  can  be 
derived  from  Darcy's  law.  The  groundwater  velocity  at  each  node  is  calculated  utilizing  an 
explicit  finite-difference  approximation  of  Darcy's  law.  The  computer  program  uses  an 
alternating-direction  implicit  procedure  to  solve  a  finite-difference  approximation  to  the 
groundwater  flow  equation,  and  it  uses  the  method  of  characteristics  (MOC)  to  solve  the 
solute  transport  equation.  MOC  uses  a  particle  tracking  procedure  to  represent  convective 
transport  and  a  two-step,  explicit  procedure  to  solve  a  finite-difference  equation  that 
describes  the  effects  of  hydrodynamic  dispersion,  fluid  sources  and  sinks,  and  divergence 
of  velocity. 

A  number  of  assumptions  are  inherent  in  the  solute  transport  model; 

1.  Darcy's  law  is  valid  and  hydraulic  head  gradients  are  the  only  significant  driving 
mechanism  for  fluid  flow. 

2.  The  porosity  and  hydraulic  conductivity  of  the  aquifer  are  constant  with  time,  and 
porosity  is  uniform  in  space. 

3.  Gradients  of  fluid  density,  viscosity,  and  temperature  do  not  affect  the  velocity 
distribution. 

4.  No  chemical  reactions  occur  that  affect  the  concentration  of  the  solute,  the  fluid 
properties,  or  the  aquifer  properties. 

5.  Ionic  and  molecular  diffusion  are  negligible  contributors  to  the  total  dispersive  flux. 

6.  Vertical  variations  in  head  and  concentrations  are  negligible. 
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7.  The  aquifer  is  homogeneous  and  isotropic  with  respect  to  the  coefficients  of 
longitudinal  and  transverse  dispersivity. 


Transport  Model  Input 

A  model  grid  of  32  columns  by  1 9  rows  with  a  250  foot  lateral  spacing  was  used. 
Specified  head  cells  were  used  at  nodes  corresponding  to  the  South  Branch  of  the  Thunder 
Bay  River,  at  nodes  along  the  eastern  boundary  of  the  model  grid  area,  and  also  at  the 
sinkhole  in  the  northeastern  portion  of  the  model.  Groundwater  elevations  measured 
during  September,  1993  were  used  as  initial  input  into  the  transport  model. 

Hydraulic  conductivity  values  were  calculated  from  slug  tests  performed  at  Alpena  CRTC 
(Engineering  Science,  1989;  Earth  Technology,  1994).  Values  of  hydraulic  conductivity 
range  from  1 2  feet/day  at  Site  4  to  278  feet/day  at  Site  3. 

Aquifer  thickness  values  'were  obtained  from  drilling  records  of  monitoring  wells  and  soil 
borings  obtained  from  the  SI  and  Rl  field  activities.  Values  listed  are  from  logs  in  which 
the  thickness  of  the  shallow  aquifer  was  clearly  discernible,  and  ranged  from  20  feet  at 
Site  5  to  65  feet  at  Site  8.  Transmissivity  values  were  calculated  by  multiplying  the 
calculated  hydraulic  conductivity  values  by  the  aquifer  thickness.  Transmissivity  ranges 
from  420  ft^/day  at  TF4-MW3  to  15,290  ft^/day  at  CG3-MW5. 

Monitoring  of  the  discharge  of  springs  into  the  sinkhole  was  performed  during  the  SI 
(Engineering  Science,  1990)  and  an  estimate  of  approximately  18,000  gallons  of  water  per 
day  discharging  into  the  sinkhole  was  calculated.  In  order  to  obtain  a  numerical  estimate 
of  discharge  into  the  sinkhole  for  the  model,  MODFLOW  (McDonald  and  Harbaugh,  1988), 
a  3-  dimensional  finite-difference  groundwater  flow  model  was  used,  MODFLOW  was 
used  because  of  its  ability  to  simulate  the  effect  of  head-dependent  groundwater  flow  into 
a  groundwater  sink  (i.e.  the  sinkhole).  This  package  was  not  available  in  MOC,  The  same 
model  parameters  and  boundary  conditions  were  used  within  MODFLOW  as  in  MOC. 

Based  upon  hydraulic  head  data  collected  in  September  1993,  discharge  from  the  shallow 
aquifer  into  the  sinkhole  is  approximately  30,000  gallons  per  day. 

The  dispersivity  of  an  aquifer  in  two  dimensions  is  described  by  the  longitudinal  dispersion, 
the  transverse  dispersion  and  the  ratio  of  the  two  (Fetter,  1993).  As  a  contaminant  plume 
moves  further  from  its  initial  location  within  the  aquifer  by  advection  with  the  groundwater 
flow,  the  plume  spreads.  The  spreading  in  the  direction  of  groundwater  flow  is  the 
longitudinal  dispersion,  the  spreading  in  the  direction  perpendicular  to  the  groundwater 
flow  is  known  as  the  transverse  dispersion  (Fetter,  1993).  The  values  of  the  dispersion 
coefficients  are  typically  derived  via  bench  scale  tests,  aquifer  tests,  or  calibration  of 
contaminant  transport  models.  Since  no  data  presently  exists  describing  dispersivity 
within  the  shallow  aquifer  beneath  the  Alpena  CRTC  and  insufficient  data  exists  to  allow 
for  derivation  of  dispersivity  via  model  calibration,  moderate  values  of  100  feet  for 
longitudinal  dispersivity  and  30  feet  for  transverse  dispersivity  were  chosen  (Gillham  and 
Cherry,  1982).  A  more  complete  description  of  the  model  is  given  in  the  report. 

Preliminary  Groundwater  Modeling  Effort,  Earth  Technology,  August  1993. 
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Model  Calibration: 


The  groundwater  flow  model  was  calibrated  with  respect  to  the  September  1993 
groundwater  elevation  measurements.  Calibration  of  the  groundwater  flow  model  was 
accomplished  by  defining  a  set  of  parameters,  boundary  conditions,  and  stresses  that 
produce  simulated  heads  and  fluxes  that  match  field-measured  values  within  a 
preestablished  range  of  error  (Anderson  and  Woessner,  1992).  In  order  to  match  field 
measured  values  for  hydraulic  head  as  determined  during  September  1993,  a  few 
modifications  were  made  to  the  preexisting  groundwater  flow  model.  These  changes 
included  updating  the  initial  head  array,  modeling  the  sinkhole  as  a  constant  head  cell  to 
account  for  the  large  gradient  changes  in  the  vicinity  of  the  sinkhole  and  including  recharge 
to  the  model  at  a  rate  of  9  inches  per  year  over  the  whole  model  area.  By  adjusting  these 
parameters,  an  acceptable  level  of  calibration  was  achieved.  An  acceptable  level  of 
calibration  was  defined  as  a  root  mean  squared  error  (RMS)  of  less  than  2  feet.  The  RMS, 
or  the  standard  deviation  is  the  average  squared  difference  in  measured  and  simulated 
heads  and  is  given  by  the  equation: 


RMS=[i/n^ 


n  =  number  of  wells 
=  measured  head 
h,  =  model  simulated  head 


Certain  portions  of  the  model  may  have  values  above  the  goal  of  2  feet  while  others  fall 
much  below  this  value.  The  RMS  represents  the  average  error  present  in  the  model. 
The  following  provides  a  summary  of  the  final  calibrated  heads  for  the  flow  model. 


Well  # 

Row,  Col 

Measured 

Head  (ft) 

Simulated 

Head  (ft) 

Measured  • 
Simulated 

Head  (ft) 

MP2MW1 

10,25 

679.69 

676.68 

3.01 

MP2MW3 

10,22 

675.51 

675.46 

0.05 

MP2MW4_5* 

11,23 

675.34 

675.70 

-0.36 

MP2MW6 

11,20 

674.86 

674.52 

0.34 

CG3MW1 

6,24 

677.38 

1 

676.69 

0.69 

CG3MW2 

7,22 

676.29 

675.78 

0.51 

CG3MW3 

9,23 

676.50 

675.98 

0.52 

CG3MW4_5* 

8,23 

676.41 

676.08 

0.33 
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Well  # 

Row,  Col 

Measured 

Head  (ft) 

Simulated 

Head  (ft) 

Measured  - 
Simulated 

Head  (ft) 

CG3MW7 

8,20 

675.64 

674.80 

0.84 

TF4MW1 

7,9 

667.23 

668.83 

-1.60 

TF4MW2 

8,10 

659.61 

664.22 

-4.61 

TF4MW3_4* 

9,10 

658.21 

660.06 

-1.85 

SF5MW1 

12,6 

674.15 

671.34 

2.81 

SF5MW2 

14,5 

675.32 

672.67 

2.65 

SF5MW3_4» 

13,6 

676.46 

671.64 

4.82 

SF5MW6 

13,5 

674.26 

672.40 

1.86 

LF6MW1 

14,9 

672.1 

670.40 

1.70 

LF6MW2 

14,8 

672.68 

670.72 

1.96 

LF6MW3 

14,10 

671.93 

670.21 

1.72 

LF6MW4 

16,7 

672.75 

671.91 

0.84 

LF6MW5 

15,7 

673.07 

671.56 

1.51 

LF6MW6 

13,10 

671.17 

669.67 

1.50 

LF6MW8 

15,9 

673.12 

670.78 

2.34 

HN8MW1 

5,22 

676.96 

675.93 

1.03 

HN8MW2 

6,19 

675.31 

674.50 

0.81 

HN8MW3_4* 

7,21 

676.01 

675.35 

0.66 

RT9MW1 

6,16 

673.06 

672.78 

0.28 

RT9MW2 

7,14 

668.21 

670.81 

-2.60 

RT9MW3 

9,15 

670.72 

671.26 

-0.54 

RT9MW4_5* 

8,14 

667.47 

670.32 

-2.85 

RT9MW6 

8,16 

670.58 

672.33 

-1.75 

S1MW2 

13,26 

677.39 

676.63 

0.76 

S1MW3 

13,27  - 

677.15 

676.98 

0.17 

S1MW1 1 

15,24 

675.72 

675.67 

0.05 
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Well  # 

Row,  Col 

Measured 

Head  (ft) 

Simulated 

Head  (ft) 

Measured  - 
Simulated 

Head  (ft) 

S1MW12 

16,25 

674.55 

675.94 

-1.39 

S1MW13 

14,24 

675.21 

675.76 

-0.55 

S1MW14 

14,25 

673.92 

676.14 

-2.22 

MP2MW2 

12,24 

675.57 

675.96 

-0.39 

*  Indicates  that  more  than  one  well  is  present  in  each  node  and  an  average  value  for 
hydraulic  head  was  used. 

Sum  of  Squared  Residuals  =  128.78/38  =  3.3891 
Root  Mean  Squared  Error  =  1 .84 

It  should  also  be  noted  that  the  model  was  calibrated  with  respect  to  the  September  1 993 
water  level  data  and  should  only  be  considered  calibrated  with  respect  to  this  data.  More 
information  on  the  water  level  fluctuation  through  time  would  be  needed  to  perform  a 
transient  calibration.  The  model  was  not  calibrated  with  respect  to  concentration  data,  but 
only  with  respect  to  hydraulic  head. 

Model  Assumptions  and  Limitations 

*  The  model  domain  consisted  only  of  the  shallow  unconfined  aquifer  (i.e.  one  layer). 

*  The  initial  head  data  input  to  the  transport  model  are  results  of  measurements  taken  in 
September  1993. 

*  Initial  concentrations  of  compounds  are  results  of  the  Round  IV  sampling  event  which 
was  conducted  from  July  to  September  1 993. 

*  Hydraulic  conductivity  values  are  the  result  of  slug  tests  performed  in  November,  1987 
and  September  1993. 

*  The  model  was  calibrated  with  respect  to  hydraulic  head  using  September  1993  water 
level  data  and  should  only  be  considered  calibrated  with  respect  to  September  1993 
water  level  data. 

*  The  flow  model  was  assumed  to  be  at  steady-state  with  respect  to  hydraulic  head. 
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Site  5 


Several  chemicals  were  detected  in  well  SF5MW1  above  MDNR  Type  A  or  Type  B  cleanup 
criteria.  These  compounds  are  1,2-  Dichloroethane,  1,4-  Dichlorobenzene,  Benzene,  and 
Styrene  at  concentrations  of  0.44,  3.5,  52,  and  1.5  ug/i,  respectively.  Chemical  data  was 
input  to  the  model  and  concentrations  of  the  contaminants  were  monitored  along  the 
eastern  shore  of  Lake  Winyah  and  at  the  sinkhole  with  respect  to  time. 
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Appendix  U:  Site  6/7  Risk  Assessment 


Table  U-1A 

a  Utilized  to  Calculate  Reasonable  Maximum  Exposure  Concentrations  for  Groundwater  -Site  6/7 

MIANG,  Alpena  CRTC,  Alpena,  Ml 
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Table  U.2A 

Data  Utilized  to  Calculate  Reasonable  Maximum  Exposure  Concentrations  for  Sediment  -  Site  6/7 

MIANG,  Alpena  CRTC,  Alpena,  Ml 
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TABLE  U-4  DAILY  IMTAKE  -  ADULT 

DERMAL  CONTACT  WITH  SEDIMENT 
FUTURE  LAND  USE  SCENARIO 


TABLE  U-6  EXPOSURE  ASSESSMENT  PARAMETERS  -  SEDIMENT  INGESTION 

Sit*  8  A  7,  Forrmr  Solid  W*«t*  Landfill  and  Firat  Fir*  Training  Araa,  Alpana  CRTC,  Atpona,  Ml 
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TABLE  U-1 1  EXPOSURE  ASSESSMENT  PARAMETERS  ■  FISH  INGESTION 

Site  6  &  7,  Former  Solid  Waste  Landfili  and  First  Fire  Training  Area,  Alpena  CRTC,  Alpena,  Ml 
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TABLE  U-1 2  DAILY  INTAKE  -  ADULT 

FISH  INGESTION 
FUTURE  LAND  USE  SCENARIO 


TABLE  U-14  DAILY  INTAKE  -  ADULT 
FISH  INGESTION 
FUTURE  LAND  USE  SCENARIO 


TABLE  U-16  EXPOSURE  ASSESSMENT  PARAMETERS  -  INGESTION  OF  SURFACE  WATER 

Site  6  &  7,  Former  Solid  Waste  Landfill  and  First  Fire  Training  Area,  Alpena  CRTC,  Alpena,  Ml 
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TABLE  U-17  DAILY  INTAKE  -  ADULT 
SURFACE  WATER  INGESTION 
FUTURE  LAND  USE  SCENARIO 


EXPOSURE  ASSESSMENT  PARAMETERS  -  DERMAL  CONTACT  WITH  THUNDERBAY  RIVER  SURFACE  WATER 
6  &  7,  Former  Solid  Waste  Landfill  and  First  Fire  Training  Area,  Alpena  CRTC,  Alpena,  Ml 


TABLE  U-22  DAILY  INTAKE  -  ADULT 

SURFACE  WATER  DERMAL  CONTACT 
FUTURE  LAND  USE  SCENARIO 
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TABLE  U-24  DAILY  INTAKE  -  ADULT 

SURFACE  WATER  DERMAL  CONTACT 
FUTURE  LAND  USE  SCENARIO 
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TABLE  U-27  SUBCHRONIC  HAZARDOUS  INDEX  ESTIMATE  ■  SEDIMENT  INGESTION 
FUTURE  LAND  USE  SCENARIO  ■  RECREATIONAL  ADULT  AND  CHILD 


FUTURE  LAND  USE  SCENARIO 


Solid  Waste  Landfill  and  First  Fire  Training  Area.  Alpena  CRTC,  Alpena,  Ml 


TABLE  U-30  CANCER  ESTIMATE  -  FISH  INGESTION 

FUTURE  LAND  USE  SCENARIO  -  ADULTS  AND  CHILDREN 


TABLE  U.31  CHRONIC  HAZARDOUS  INDEX  ESTIMATE  -  FISH  INGESTION 
FUTURE  LAND  USE  SCENARIO  -  ADULTS  AND  CHILDREN 


TABLE  U-32  CANCER  ESTIMATE  -  SURFACE  WATER  INGESTION 
FUTURE  LAND  USE  SCENARIO  -  ADULTS  AND  CHILDREN 


TABLE  U-33  CHRONIC  HAZARDOUS  INDEX  ESTIMATE  -  SURFACE  WATER  INGESTION 
FUTURE  LAND  USE  SCENARIO  -  ADULTS  AND  CHILDREN 


TABLE  U-34  CANCER  ESTJMATE  •  SURFACE  WATER  DERMAL 
FUTURE  LAND  USE  SCENARIO  -  ADULTS  AND  CHILDREN 


TABLE  U-35  CHRONIC  HAZARDOUS  INDEX  ESTIMATE  -  SURFACE  WATER  DERMAL 
FUTURE  LAND  USE  SCENARIO  -  ADULTS  AND  CHILDREN 


Groundwater  Solute  Transport  Model  Data  -  Site  6/7 
Theoretical  Background 


A  two-dimensional  Method  of  Characteristics  (MOC)  solute  transport  model  (Konikow  and 
Bredehoeft,  1989)  was  used  for  preliminary  examination  of  contaminant  migration  within 
the  shallow  aquifer  beneath  the  Alpena  CRTC.  The  model  is  designed  to  calculate 
transient  changes  in  solute  concentrations  within  groundwater  by  simultaneously  solving 
partial  differential  equations  describing  groundwater  flow  and  transport  and  computes  the 
change  in  a  chemicals  concentration  over  time.  Changes  in  chemical  concentrations  over 
time  are  caused  by  the  processes  of  convective  transport,  hydrodynamic  dispersion,  and 
mixing  from  fluid  sources.  This  model  couples  the  groundwater  flow  equation  with  solute 
transport  equations. 

The  flow  equation  can  be  approximated  by  an  implicit  finite-difference  equation.  The 
model  area  is  discretized  into  a  rectangular  grid  with  each  square  being  a  node.  The  finite 
difference  equation  is  solved  numerically  for  each  node  in  the  grid  using  an  iterative 
alternating-direction  implicit  (ADI)  procedure. 

After  the  hydraulic  head  distribution  is  calculated,  the  velocity  of  groundwater  flow  can  be 
computed  at  each  node.  The  expression  for  average  velocity  of  groundwater  can  be 
derived  from  Darcy's  law.  The  groundwater  velocity  at  each  node  is  calculated  utilizing  an 
explicit  finite-difference  approximation  of  Darcy's  law.  The  computer  program  uses  an 
alternating-direction  implicit  procedure  to  solve  a  finite-difference  approximation  to  the 
groundwater  flow  equation,  and  it  uses  the  method  of  characteristics  (MOC)  to  solve  the 
solute  transport  equation.  MOC  uses  a  particle  tracking  procedure  to  represent  convective 
transport  and  a  two-step,  explicit  procedure  to  solve  a  finite-difference  equation  that 
describes  the  effects  of  hydrodynamic  dispersion,  fluid  sources  and  sinks,  and  divergence 
of  velocity. 

A  number  of  assumptions  are  inherent  in  the  solute  transport  model: 

1 .  Darcy's  law  is  valid  and  hydraulic  head  gradients  are  the  only  significant  driving 
mechanism  for  fluid  flow. 

2.  The  porosity  and  hydraulic  conductivity  of  the  aquifer  are  constant  with  time,  and 
porosity  is  uniform  in  space. 

3.  Gradients  of  fluid  density,  viscosity,  and  temperature  do  not  affect  the  velocity 
distribution. 

4.  No  chemical  reactions  occur  that  affect  the  concentration  of  the  solute,  the  fluid 
properties,  or  the  aquifer  4)roperties. 

5.  Ionic  and  molecular  diffusion  are  negligible  contributors  to  the  total  dispersive  flux. 

6.  Vertical  variations  in  head  and  concentrations  are  negligible. 

Site  6  /  gwmod.sS  U-29 


7.  The  aquifer  is  homogeneous  and  isotropic  with  respect  to  the  coefficients  of 
longitudinal  and  transverse  dispersivity. 


Transport  Model  Input 

A  model  grid  of  32  columns  by  19  rows  with  a  250  foot  lateral  spacing  was  used. 
Specified  head  cells  were  used  at  nodes  corresponding  to  the  South  Branch  of  the  Thunder 
Bay  River,  at  nodes  along  the  eastern  boundary  of  the  model  grid  area,  and  also  at  the 
sinkhole  in  the  northeastern  portion  of  the  model.  Groundwater  elevations  measured 
during  September,  1993  were  used  as  initial  input  into  the  transport  model. 

Hydraulic  conductivity  values  were  calculated  from  slug  tests  performed  at  Alpena  CRTC 
(Engineering  Science,  1989;  Earth  Technology,  1994).  Values  of  hydraulic  conductivity 
range  from  1 2  feet/day  at  Site  4  to  278  feet/day  at  Site  3. 

Aquifer  thickness  values  were  obtained  from  drilling  records  of  monitoring  wells  and  soil 
borings  obtained  from  the  SI  and  Rl  field  activities.  Values  listed  are  from  logs  in  which 
the  thickness  of  the  shallow  aquifer  was  clearly  discernible,  and  ranged  from  20  feet  at 
Site  5  to  65  feet  at  Site  8.  Transmissivity  values  were  calculated  by  multiplying  the 
calculated  hydraulic  conductivity  values  by  the  aquifer  thickness.  Transmissivity  ranges 
from  420  ftVday  at  TF4-MW3  to  15,290  ft^/day  at  CG3-MW5. 

Monitoring  of  the  discharge  of  springs  into  the  sinkhole  was  performed  during  the  SI 
(Engineering  Science,  1990)  and  an  estimate  of  approximately  18,000  gallons  of  water  per 
day  discharging  into  the  sinkhole  was  calculated.  In  order  to  obtain  a  numerical  estimate 
of  discharge  into  the  sinkhole  for  the  model,  MODFLOW  (McDonald  and  Harbaugh,  1988), 
a  3-  dimensional  finite-difference  groundwater  flow  model  was  used.  MODFLOW  was 
used  because  of  its  ability  to  simulate  the  effect  of  head-dependent  groundwater  flow  into 
a  groundwater  sink  (i.e.  the  sinkhole).  This  package  was  not  available  in  MOC.  The  same 
model  parameters  and  boundary  conditions  were  used  within  MODFLOW  as  in  MOC. 

Based  upon  hydraulic  head  data  collected  in  September  1993,  discharge  from  the  shallow 
aquifer  into  the  sinkhole  is  approximately  30,000  gallons  per  day. 

The  dispersivity  of  an  aquifer  in  two  dimensions  is  described  by  the  longitudinal  dispersion, 
the  transverse  dispersion  and  the  ratio  of  the  two  (Fetter,  1993).  As  a  contaminant  plume 
moves  further  from  its  initial  location  within  the  aquifer  by  advection  with  the  groundwater 
flow,  the  plume  spreads.  The  spreading  in  the  direction  of  groundwater  flow  is  the 
longitudinal  dispersion,  the  spreading  in  the  direction  perpendicular  to  the  groundwater 
flow  is  known  as  the  transverse  dispersion  (Fetter,  1993).  The  values  of  the  dispersion 
coefficients  are  typically  derived  via  bench  scale  tests,  aquifer  tests,  or  calibration  of 
contaminant  transport  models.  Since  no  data  presently  exists  describing  dispersivity 
within  the  shallow  aquifer  beneath  the  Alpena  CRTC  and  insufficient  data  exists  to  allow 
for  derivation  of  dispersivity  via  model  calibration,  moderate  values  of  100  feet  for 
longitudinal  dispersivity  and  30  feet  for  transverse  dispersivity  were  chosen  (Gillham  and 
Cherry,  1982).  A  more  complete  description  of  the  model  is  given  in  the  report. 

Preliminary  Groundwater  Modeling  Effort,  Earth  Technology,  August  1993. 
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Model  Calibration: 


The  groundwater  flow  model  was  calibrated  with  respect  to  the  September  1993 
groundwater  elevation  measurements.  Calibration  of  the  groundwater  flow  model  was 
accomplished  by  defining  a  set  of  parameters,  boundary  conditions,  and  stresses  that 
produce  simulated  heads  and  fluxes  that  match  field-measured  values  within  a 
preestablished  range  of  error  (Anderson  and  Woessner,  1992).  In  order  to  match  field 
measured  values  for  hydraulic  head  as  determined  during  September  1993,  a  few 
modifications  were  made  to  the  preexisting  groundwater  flow  model.  These  changes 
included  updating  the  initial  head  array,  modeling  the  sinkhole  as  a  constant  head  cell  to 
account  for  the  large  gradient  changes  in  the  vicinity  of  the  sinkhole  and  including  recharge 
to  the  model  at  a  rate  of  9  inches  per  year  over  the  whole  model  area.  By  adjusting  these 
parameters,  an  acceptable  level  of  calibration  was  achieved.  An  acceptable  level  of 
calibration  was  defined  as  a  root  mean  squared  error  (RMS)  of  less  than  2  feet.  The  RMS, 
or  the  standard  deviation  is  the  average  squared  difference  in  measured  and  simulated 
heads  and  is  given  by  the  equation: 


n  =  number  of  wells 
h„  =  measured  head 
h,  =  model  simulated  head 

Certain  portions  of  the  model  may  have  values  above  the  goal  of  2  feet  while  others  fall 
much  below  this  value.  The  RMS  represents  the  average  error  present  in  the  model. 

The  following  provides  a  summary  of  the  final  calibrated  heads  for  the  flow  model. 


Well  # 

Row,  Col 

Measured 

Head  (ft) 

Simulated 

Head  (ft) 

Measured  - 
Simulated 

Head  (ft) 

CG3MW7 

8,20 

675.64 

674.80 

0.84 

TF4MW1 

7,9 

667.23 

668.83 

-1.60 

TF4MW2 

8,10 

659.61 

664.22 

-4.61 

TF4MW3_4* 

9,10 

658.21 

660.06 

-1.85 

SF5MW1 

12,6 

674.15 

671.34 

2.81 

SF5MW2 

14,5 

675.32 

672.67 

2.65 

SF5MW3_4* 

13,6 

676.46 

671.64 

4.82 

SF5MW6 

13,5 

674.26 

672.40 

1.86 

LF6MW1 

14,9 

672.1 

670.40 

1.70 

LF6MW2 

14,8 

672.68 

670.72 

1.96 

LF6MW3 

14,10 

671.93 

670.21 

1.72 

LF6MW4 

16,7 

672.75 

671.91 

0.84 

LF6MW5 

15,7 

673.07 

671.56 

1.51 

LF6MW6 

13,10 

671.17 

669.67 

1.50 

LF6MW8 

15,9 

673.12 

670.78 

2.34 

HN8MW1 

5,22 

676.96 

675.93 

1.03 

HN8MW2 

6,19 

675.31 

674.50 

0.81 

HN8MW3_4* 

7,21 

676.01 

675.35 

0.66 

RT9MW1 

6,16 

673.06 

672.78 

0.28 

RT9MW2 

7,14 

668.21 

670.81 

-2.60 

RT9MW3 

9,15 

670.72 

671.26 

-0.54 

RT9MW4_5* 

8,14 

667.47 

670.32 

-2.85 

RT9MW6 

8,16 

670.58 

672.33 

-1.75 

S1MW2 

13,26 

677.39 

676.63 

0.76 

S1MW3 

13,27  " 

677.15 

676.98 

0.17 

S1MW1  1 

15,24 

675.72 

675.67 

0.05 

S1MW12 

16,25 

674.55 

675.94 

-1.39 
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Well  # 

Row,  Col 

Measured 

Head  (ft) 

Simulated 

Head  (ft) 

Measured  - 
Simulated 

Head  (ft) 

S1MW13 

14,24 

675.21 

675.76 

-0.55 

S1MW14 

14,25 

673.92 

676.14 

-2.22 

MP2MW2 

12,24 

675.57 

675.96 

-0.39 

*  Indicates  that  more  than  one  well  is  present  in  each  node  and  an  average  value  for 
hydraulic  head  was  used. 

Sum  of  Squared  Residuals  =  128.78/38  =  3.3891 
Root  Mean  Squared  Error  =  1.84 


It  should  also  be  noted  that  the  model  was  calibrated  with  respect  to  the  September  1993 
water  level  data  and  should  only  be  considered  calibrated  with  respect  to  this  data.  More 
information  on  the  water  level  fluctuation  through  time  would  be  needed  to  perform  a 
transient  calibration.  The  model  was  not  calibrated  with  respect  to  concentration  data,  but 
only  with  respect  to  hydraulic  head. 

Model  Assumptions  and  Limitations 

*  The  model  domain  consisted  only  of  the  shallow  unconfined  aquifer  (i.e.  one  layer). 

*  The  initial  head  data  input  to  the  transport  model  are  results  of  measurements  taken  in 
September  1993. 

*  Initial  concentrations  of  compounds  are  results  of  the  Round  IV  sampling  event  which 
was  conducted  from  July  to  September  1993. 

*  Hydraulic  conductivity  values  are  the  result  of  slug  tests  performed  in  November,  1987 
and  September  1993. 

*  The  model  was  calibrated  with  respect  to  hydraulic  head  using  September  1 993  water 
level  data  and  should  only  be  considered  calibrated  with  respect  to  September  1993 
water  level  data. 

*  The  flow  model  was  assumed  to  be  at  steady-state  with  respect  to  hydraulic  head. 
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Site  6 


Carbon  Tetrachloride  was  present  in  LF6MW3  above  MDNR  Type  B  cleanup  criteria  at  a 
concentration  of  1 .2  ug/1.  This  data  was  input  to  the  model  and  concentrations  were 
monitored  along  the  backwater  area  of  Lake  Winyah  and  at  the  sinkhole  with  respect  to 
time. 
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Appendix  V:  Site  8  Risk  Assessment 


Table  V-1A 

Data  Utilized  to  Calculate  Reasonable  Maximum  Exposure  Concentrations  for  Groundwater  -  Site  8 

MIANG,  Alpena  CRTC,  Alpena,  Ml 


SITE8V.WB1/ 940126 


MIANG,  Alpena  CRTC,  Alpena,  Ml 


SITE8V.WB1/940126 


TABLE  V-1  EXPOSURE  ASSESSMENT  PARAMETERS  -  INGESTION  OF  GROUNDWATER 
Site  8,  Former  Site  of  Hangar  9,  Alpena  CRTC,  Alpena,  Ml 


1/940127 


TABLE  V-2  DAILY  INTAKE  -  ADULT 

INGESTION  OF  GROUNDWATER 
FUTURE  LAND  USE  SCENARIO 


V-4 


Carcinogenic  Effect  Exposure  Assessment 

Site  8,  Former  Site  of  Hangar  9,  Alpena  CRTC,  Alpena,  Ml 


SITE8A.WB1/940127 


TABLE  V-6  EXPOSURE  ASSESSMENT  PARAMETERS -INHALATION  OF  GROUNDWATER 
Site  8,  Former  Site  of  Hangar  9,  Alpena  CRTC,  Alpena,  Ml 


c  c 

mS  C 
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o  o  >  2  H  c 
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TABLE  V-7  DAILY  INTAKE  .  ADULT 

INHALATION  OF  GROUNDWATER 
FUTURE  LAND  USE  SCENARIO 


SITE8A.WB1/940127 


TABLE  V-9  DAILY  INTAKE  -  ADULT 

INHALATION  OF  GROUNDWATER 
FUTURE  LAND  USE  SCENARIO 


V-8 


TABLE  V<11  GROUNDWATER  INHALATION  MODEL  CALCULATIONS 
Site  8,  Former  Site  of  Hangar  9,  Alpena  CRTC,  Alpena,  Ml 


SITE8A.WB1/ 940127 


SITE8A.WB1/ 940127 


TABLE  V-15  DAILY  INTAKE  -  ADULT 

GROUNDWATER  DERMAL  CONTACT 
FUTURE  LAND  USE  SCENARIO 


V-12 


TABLE  V-18  DAILY  INTAKE -ADULT 

DERMAL  CONTACT  WITH  SOILS 

CURRENT  LAND  USE  SCENARIO  -  ON-SITE  EMPLOYEE 


TABLE  V.19  DAILY  INTAKE  •  ADULT 

DERMAL  CONTACT  WITH  SOILS 

CURRENT  LAND  USE  SCENARIO  .  ON-SITE  EMPLOYEE 


SITE8AWB1/ 940127 


TABLE  V-20  DAILY  INTAKE  -  ADULT 

DERMAL  CONTACT  WITH  SOILS 

FUTURE  LAND  USE  SCENARIO  -  EXCAVATION  WORKERS 


S.  «  ° 
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TABLE  V-21  DAILY  INTAKE -ADULT 

DERMAL  CONTACT  WITH  SOILS 

FUTURE  LAND  USE  SCENARIO  -  EXCAVATION  WORKERS 


TABLE  V.22  DAILY  INTAKE  -  ADULTS 
DERMAL  CONTACT  WITH  SOILS 

FUTURE  LAND  USE  SCENARIO  -  RECREATIONAL  ADULT 


V-18 


940122 


TABLE  V.24  DAILY  tNTAKE  >  ADULTS 
DERMAL  CONTACT  WITH  SOILS 

FUTURE  LAND  USE  SCENARIO  -  RECREATIONAL  ADULT 


SITE8A.WB1/ 940127 


EXPOSURE  ASSESSMENT  PARAMETERS  -  SOIL  INGESTION 
i,  Fortmr  Site  of  Hangar  9,  Alpena  CRTC.  Alpena.  MI 
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Noncarcinogenic  Effect  Exposure  Assessment 

Site  8,  Former  Site  of  Hangar  9.  Atpena  CRTC,  Alpena,  Ml 
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TABLE  V-30  DAILY  INTAKE  -  ADULT 
SOIL  INGESTION 


I 


V-25 
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Carcinogenic  Effect  Exposure  Assessment 

Site  8,  Former  Site  of  Hangar  9,  Alpena  CRTC,  Alpena,  Ml 


SITE8A.WB1/ 940127 


Noncarcinogenic  Effect  Exposure  Assessment 

Site  8,  Former  Site  of  Hangar  9,  Alpena  CRTC,  Alpena  Ml 


SITE8A.WB1/ 940127 


TABLE  V-37  DAILY  INTAKE  -  ADULT 
SOIL  INHALATION 

FUTURE  LAND  USE  SCENARIO  -  EXCAVATION  WORKERS 


TABLE  V-38  SOIL  INHALATION  MODEL  CALCULATIONS 

Site  8,  Former  Site  of  Hangar  9,  Alpena  CRTC,  Alpena,  Ml 


SITE8A.WB  1/9401 27 


TABLE  V-39  SOIL  INHALTION  MODEL  CALCULATIONS 

Site  8,  Former  Site  of  Hangar  9,  Alpena  CRTC,  Alpena,  Ml 


TABLE  V^O  CANCER  ESTtMATE -GROUNDWATER  INGESTION 
FUTURE  LAND  USE  SCENARIO  -  ADULTS  AND  CHILDREN 


SITESA.WBI/  940127 


TABLE  V-4a  CANCER  ESTIMATE  -  GROUNDWATER  INHALATION 
FUTURE  LAND  USE  SCENARIO  -  ADULTS  AND  CHILDREN 


TABLE  V-44  CHRONIC  HAZARDOUS  INDEX  ESTIMATE  *  GROUNDWATER  INHALATION 
FUTURE  LAND  USE  SCENARIO  -  ADULTS  AND  CHILDREN 


SITE  8A.WB  1/940127 


TABLE  V-45  CANCER  ESTIMATE  -  DERMAL  CONTACT  WITH  GROUNDWATER 
FUTURE  LAND  USE  SCENARIO  -  ADULTS  AND  CHILDREN 


V-36 


V-4G  CHRONIC  HAZARDOUS  INDEX  ESTIMATE  -  DERMAL  CONTACT  WITH  GROUNDWATER 
FUTURE  LAND  USE  SCENARIO  -  ADULTS  AND  CHILDREN 


SITE8A.WB1/ 940127 


TABLE  V-»7  CANCER  ESTIMATE-  DERMAL  CONTACT  WITH  SCrtLS 


TABLe  V-4«  CANCBR  BSnMATB-  DERMAL  CONTACT  WITH  SOILS  -  OM.«TE  EMPLOYEE 


TAiLi  V.M  SUBCHRONIC  HAZARt>OUS  INDEX  ESTIMATE  -  DERMAL  CONTACT  WITH  SOILS  -  ON-StTE  EMELOVEE 


TABLE  V-55  CANCER  ESTIMATE  -  SOIL  INHALATION 

FUTURE  LAND  USE  SCENARIO  .  EXCAVATION  WORKER 


SITEfiAWBI/ 940127 


TABLE  V-56  SUBCHRONIC  HAZARDOUS  INDEX  ESTIMATE  -  SOIL  INHALATION 
FUTURE  LAND  USE  SCENARIO  -  EXCAVATION  WORKER 


Groundwater  Solute  Transport  Model  Data  -  Site  8 
Theoretical  Background 

A  two-dimensional  Method  of  Characteristics  (MOC)  solute  transport  model  (Konikow  and 
Bredehoeft,  1 989)  was  used  for  preliminary  examination  of  contaminant  migration  within 
the  shallow  aquifer  beneath  the  Alpena  CRTC.  The  model  is  designed  to  calculate 
transient  changes  in  solute  concentrations  within  groundwater  by  simultaneously  solving 
partial  differential  equations  describing  groundwater  flow  and  transport  and  computes  the 
change  in  a  chemicals  concentration  over  time.  Changes  in  chemical  concentrations  over 
time  are  caused  by  the  processes  of  convective  transport,  hydrodynamic  dispersion,  and 
mixing  from  fluid  sources.  This  model  couples  the  groundwater  flow  equation  with  solute 
transport  equations. 

The  flow  equation  can  be  approximated  by  an  implicit  finite-difference  equation.  The 
model  area  is  discretized  into  a  rectangular  grid  with  each  square  being  a  node.  The  finite 
difference  equation  is  solved  numerically  for  each  node  in  the  grid  using  an  iterative 
alternating-direction  implicit  (ADI)  procedure. 

After  the  hydraulic  head  distribution  is  calculated,  the  velocity  of  groundwater  flow  can  be 
computed  at  each  node.  The  expression  for  average  velocity  of  groundwater  can  be 
derived  from  Darcy's  law.  The  groundwater  velocity  at  each  node  is  calculated  utilizing  an 
explicit  finite-difference  approximation  of  Darcy's  law.  The  computer  program  uses  an 
alternating-direction  implicit  procedure  to  solve  a  finite-difference  approximation  to  the 
groundwater  flow  equation,  and  it  uses  the  method  of  characteristics  (MOC)  to  solve  the 
solute  transport  equation.  MOC  uses  a  particle  tracking  procedure  to  represent  convective 
transport  and  a  two-step,  explicit  procedure  to  solve  a  finite-difference  equation  that 
describes  the  effects  of  hydrodynamic  dispersion,  fluid  sources  and  sinks,  and  divergence 
of  velocity. 

A  number  of  assumptions  are  inherent  in  the  solute  transport  model: 

1 .  Darcy's  law  is  valid  and  hydraulic  head  gradients  are  the  only  significant  driving 
mechanism  for  fluid  flow. 

2.  The  porosity  and  hydraulic  conductivity  of  the  aquifer  are  constant  with  time,  and 
porosity  is  uniform  in  space. 

3.  Gradients  of  fluid  density,  viscosity,  and  temperature  do  not  affect  the  velocity 
distribution. 

4.  No  chemical  reactions  occur  that  affect  the  concentration  of  the  solute,  the  fluid 
properties,  or  the  aquifer  properties. 

5.  Ionic  and  molecular  diffusion  are  negligible  contributors  to  the  total  dispersive  flux. 

6.  Vertical  variations  in  head  and  concentrations  are  negligible. 


Site  8  /  gwmod.sS 


V-48 


7.  The  aquifer  is  homogeneous  and  isotropic  with  respect  to  the  coefficients  of 
longitudinal  and  transverse  dispersivity. 


Transport  Model  Input 

A  model  grid  of  32  columns  by  19  rows  with  a  250  foot  lateral  spacing  was  used. 
Specified  head  cells  were  used  at  nodes  corresponding  to  the  South  Branch  of  the  Thunder 
Bay  River,  at  nodes  along  the  eastern  boundary  of  the  model  grid  area,  and  also  at  the 
sinkhole  in  the  northeastern  portion  of  the  model.  Groundwater  elevations  measured 
during  September,  1 993  were  used  as  initial  input  into  the  transport  model. 

Hydraulic  conductivity  values  were  calculated  from  slug  tests  performed  at  Alpena  CRTC 
(Engineering  Science,  1989;  Earth  Technology,  1994).  Values  of  hydraulic  conductivity 
range  from  12  feet/day  at  Site  4  to  278  feet/day  at  Site  3. 

Aquifer  thickness  values  were  obtained  from  drilling  records  of  monitoring  wells  and  soil 
borings  obtained  from  the  SI  and  Rl  field  activities.  Values  listed  are  from  logs  in  which 
the  thickness  of  the  shallow  aquifer  was  clearly  discernible,  and  ranged  from  20  feet  at 
Site  5  to  65  feet  at  Site  8.  Transmissivity  values  were  calculated  by  multiplying  the 
calculated  hydraulic  conductivity  values  by  the  aquifer  thickness.  Transmissivity  ranges 
from  420  ftVday  at  TF4-MW3  to  1  5,290  ft^/day  at  CG3-MW5. 

Monitoring  of  the  discharge  of  springs  into  the  sinkhole  was  performed  during  the  SI 
(Engineering  Science,  1990)  and  an  estimate  of  approximately  18,000  gallons  of  water  per 
day  discharging  into  the  sinkhole  was  calculated.  In  order  to  obtain  a  numerical  estimate 
of  discharge  into  the  sinkhole  for  the  model,  MODFLOW  (McDonald  and  Harbaugh,  1988), 
a  3-  dimensional  finite-difference  groundwater  flow  model  was  used.  MODFLOW  was 
used  because  of  its  ability  to  simulate  the  effect  of  head-dependent  groundwater  flow  into 
a  groundwater  sink  (i.e.  the  sinkhole).  This  package  was  not  available  in  MOC.  The  same 
model  parameters  and  boundary  conditions  were  used  within  MODFLOW  as  in  MOC. 

Based  upon  hydraulic  head  data  collected  in  September  1993,  discharge  from  the  shallow 
aquifer  into  the  sinkhole  is  approximately  30,000  gallons  per  day. 

The  dispersivity  of  an  aquifer  in  two  dimensions  is  described  by  the  longitudinal  dispersion, 
the  transverse  dispersion  and  the  ratio  of  the  two  (Fetter,  1993).  As  a  contaminant  plume 
moves  further  from  its  initial  location  within  the  aquifer  by  advection  with  the  groundwater 
flow,  the  plume  spreads.  The  spreading  in  the  direction  of  groundwater  flow  is  the 
longitudinal  dispersion,  the  spreading  in  the  direction  perpendicular  to  the  groundwater 
flow  is  known  as  the  transverse  dispersion  (Fetter,  1993).  The  values  of  the  dispersion 
coefficients  are  typically  derived  via  bench  scale  tests,  aquifer  tests,  or  calibration  of 
contaminant  transport  models.  Since  no  data  presently  exists  describing  dispersivity 
within  the  shallow  aquifer  beneath  the  Alpena  CRTC  and  insufficient  data  exists  to  allow 
for  derivation  of  dispersivity  via  model  calibration,  moderate  values  of  100  feet  for 
longitudinal  dispersivity  and  30  feet  for  transverse  dispersivity  were  chosen  (Gillham  and 
Cherry,  1982).  A  more  complete  description  of  the  model  is  given  in  the  report. 

Preliminary  Groundwater  Modeling  Effort,  Earth  Technology,  August  1993. 
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Model  Calibration: 


The  groundwater  flow  model  was  calibrated  with  respect  to  the  September  1993 
groundwater  elevation  measurements.  Calibration  of  the  groundwater  flow  model  was 
accomplished  by  defining  a  set  of  parameters,  boundary  conditions,  and  stresses  that 
produce  simulated  heads  and  fluxes  that  match  field-measured  values  within  a 
preestablished  range  of  error  (Anderson  and  Woessner,  1992).  In  order  to  match  field 
measured  values  for  hydraulic  head  as  determined  during  September  1993,  a  few 
modifications  were  made  to  the  preexisting  groundwater  flow  model.  These  changes 
included  updating  the  initial  head  array,  modeling  the  sinkhole  as  a  constant  head  cell  to 
account  for  the  large  gradient  changes  in  the  vicinity  of  the  sinkhole  and  including  recharge 
to  the  model  at  a  rate  of  9  inches  per  year  over  the  whole  model  area.  By  adjusting  these 
parameters,  an  acceptable  level  of  calibration  was  achieved.  An  acceptable  level  of 
calibration  was  defined  as  a  root  mean  squared  error  (RMS)  of  less  than  2  feet.  The  RMS, 
or  the  standard  deviation  is  the  average  squared  difference  in  measured  and  simulated 
heads  and  is  given  by  the  equation: 


n  =  number  of  wells 
h,^  =  measured  head 
h,  =  model  simulated  head 


Certain  portions  of  the  model  may  have  values  above  the  goal  of  2  feet  while  others  fall 
much  below  this  value.  The  RMS  represents  the  average  error  present  in  the  model. 

The  following  provides  a  summary  of  the  final  calibrated  heads  for  the  flow  model. 


Well  # 

Row,  Col 

Measured 

Head  (ft) 

Simulated 

Head  (ft) 

Measured  - 
Simulated 

Head  (ft) 

MP2MW1 

10,25 

679.69 

676.68 

3.01 

MP2MW3 

10,22 

675.51 

675.46 

0.05 

MP2MW4_5* 

11,23 

675.34 

675.70 

-0.36 

MP2MW6 

11,20 

674.86 

674.52 

0.34 

CG3MW1 

6,24 

677.38 

^  676.69 

0.69 

CG3MW2 

7,22 

676.29 

675.78 

0.51 

CG3MW3 

9,23 

676.50 

675.98 

0.52 

CG3MW4_5* 

8,23 

676.41 

676.08 

0.33 
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Well  # 

Row,  Col 

Measured 

Head  (ft) 

Simulated 

Head  (ft) 

Measured  - 
Simulated 

Head  (ft) 

CG3MW7 

8,20 

675.64 

674.80 

0.84 

TF4MW1 

7,9 

667.23 

668.83 

-1.60 

TF4MW2 

8,10 

659.61 

664.22 

-4.61 

TF4MW3_4*' 

9,10 

658.21 

660.06 

-1.85 

SF5MW1 

12,6 

674.15 

671.34 

2.81 

SF5MW2 

14,5 

675.32 

672.67 

2.65 

SF5MW3_4* 

13,6 

676.46 

671.64 

4.82 

SF5MW6 

13,5 

674.26 

672.40 

1.86 

LF6MW1 

14,9 

672.1 

670.40 

1.70 

LF6MW2 

14,8 

672.68 

670.72 

1.96 

LF6MW3 

14,10 

671.93 

670.21 

1.72 

LF6MW4 

16,7 

672.75 

671.91 

0.84 

LF6MW5 

15,7 

673.07 

671.56 

1.51 

LF6MW6 

13,10 

671.17 

669.67 

1.50 

LF6MW8 

15,9 

673.12 

670.78 

2.34 

HN8MW1 

5,22 

676.96 

675.93 

1.03 

HN8MW2 

6,19 

675.31 

674.50 

0.81 

HN8MW3_4* 

7,21 

676.01 

675.35 

0.66 

RTS  MW  1 

6,16 

673.06 

672.78 

0.28 

RT9MW2 

7,14 

668.21 

670.81 

-2.60 

RT9MW3 

9,15 

670.72 

671.26 

-0.54 

RT9MW4_5* 

8,14 

667.47 

670.32 

-2.85 

RT9MW6 

8,16 

670.58 

672.33 

-1.75 

S1MW2 

13,26 

677.39 

676.63 

0.76 

S1MW3 

13,27 

677.15 

676.98 

0.17 

S1MW11 

15,24 

675.72 

675.67 

0.05 

S1MW12 

16,25 

674.55 

675.94 

-1.39 

I 
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Well  # 

Row,  Col 

Measured 

Head  (ft) 

Simulated 

Head  (ft) 

Measured  - 
Simulated 

Head  (ft) 

S1MW13 

14,24 

675.21 

675.76 

-0.55 

S1MW14 

14,25 

673.92 

676.14 

-2.22 

MP2MW2 

12,24 

675.57 

675.96 

-0.39 

*  Indicates  that  mors  than  one  well  is  present  in  each  node  and  an  av6rage  value  for 
hydraulic  head  was  used. 

Sum  of  Squared  Residuals  =  128.78/38  =  3.3891 
Root  Mean  Squared  Error  =  1.84 


It  should  also  be  noted  that  the  model  was  calibrated  with  respect  to  the  September  1993 
water  level  data  and  should  only  be  considered  calibrated  with  respect  to  this  data.  More 
information  on  the  water  level  fluctuation  through  time  would  be  needed  to  perform  a 
transient  calibration.  The  model  was  not  calibrated  with  respect  to  concentration  data,  but 
only  with  respect  to  hydraulic  head. 

Model  Assumptions  and  Limitations 

*  The  model  domain  consisted  only  of  the  shallow  unconfined  aquifer  (i.e.  one  layer). 

*  The  initial  head  data  input  to  the  transport  model  are  results  of  measurements  taken  in 
September  1 993. 

*  Initial  concentrations  of  compounds  are  results  of  the  Round  IV  sampling  event  which 
was  conducted  from  July  to  September  1993. 

*  Hydraulic  conductivity  values  are  the  result  of  slug  tests  performed  in  November.  1987 
and  September  1993. 

*  The  model  was  calibrated  with  respect  to  hydraulic  head  using  September  1993  water 
level  data  and  should  only  be  considered  calibrated  with  respect  to  September  1993 
water  level  data. 

*  The  flow  model  was  assumed  to  be  at  steady-state  with  respect  to  hydraulic  head. 
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Site  8 


PCE  was  present  in  well  HN8MW8  above  MDNR  Type  B  cleanup  criteria  at  a  concentration 
of  1 .2  ug/l.  This  data  was  input  to  the  model  and  concentrations  at  the  sinkhole  were 
monitored  with  respect  to  time. 
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Appendix  W:  Site  9  Risk  Assessment 
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TABLE  W-1  EXPOSURE  ASSESSMENT  PARAMETERS  -  INGESTION  OF  GROUNDWATER 
Site  9,  Radar  Tower,  AJpena  CRTC,  Alpena,  Ml 
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TABLE  W-2  DAILY  INTAKE  -  ADULT 

INGESTION  OF  GROUNDWATER 
FUTURE  LAND  USE  SCENARIO 
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TABLE  W-4  DAILY  INTAKE  -  ADULT 

INGESTION  OF  GROUNDWATER 
FUTURE  LAND  USE  SCENARIO 
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TABLE  W-6  CALCULATIONS  FOR  GROUNDWATER  INHALATION 
Site  9,  Radar  Tower,  Alpena  CRTC,  Alpena,  Ml 
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TABLE  W-10  DAILY  INTAKE  -  ADULT 

INHALATION  OF  GROUNDWATER 
FUTURE  LAND  USE  SCENARIO 


TABLE  W-12  EXPOSURE  ASSESSMENT  PARAMETERS  -  DERMAL  CONTACT  WITH  GROUNDWATER 
Site  9,  Radar  Tower,  Alpena  CRTC,  Alpena,  Ml 


TABLE  W  13  DAILY  INTAKE  -  ADULT 
GROUNDWATER  DERMAL  CO^ 
FUTURE  LAND  USE  SCENARIO 
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TABLE  W-15  DAILY  INTAKE  -  ADULT 
GROUNDWATER  DERMAL  COM 
FUTURE  LAND  USE  SCENARIO 
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TABLE  W-17  CANCER  ESTIMATE  -  GROUNDWATER  INGESTION 
FUTURE  LAND  USE  SCENARIO  -  ADULTS  AND  CHILDREN 


TABLE  W  ia  CHRONIC  HAZARDOUS  INDEX  ESTIMATE  -  GROUNDWATER  INGESTION 
FUTURE  LAND  USE  SCENARIO  -  ADULTS  AND  CHILDREN 


TABLE  W-20CHRONrC  HAZARDOUS  INDEX  ESTIMATE  -  GROUNDWATER  INHALATION 
FUTURE  LAND  USE  SCENARIO  -  ADULTS  AND  CHILDREN 


TABLE  W-21  CANCER  ESTIMATE  -  DERMAL  CONTACT  WITH  GROUNDWATER 
FUTURE  LAND  USE  SCENARIO  •  ADULTS  AND  CHILDREN 


TABLE  W  22  CHRONIC  HAZARDOUS  INDEX  ESTIMATE  -  DERMAL  CONTACT  WITH  GROUNDWATER 


Groundwater  Solute  Transport  Model  Data  -  Site  9 
Theoretical  Background 

A  two-dimensional  Method  of  Characteristics  (MOC)  solute  transport  model  (Konikow  and 
Bredehoeft,  1989)  was  used  for  preliminary  examination  of  contaminant  migration  within 
the  shallow  aquifer  beneath  the  Alpena  CRTC.  The  model  is  designed  to  calculate 
transient  changes  in  solute  concentrations  within  groundwater  by  simultaneously  solving 
partial  differential  equations  describing  groundwater  flow  and  transport  and  computes  the 
change  in  a  chemicals  concentration  over  time.  Changes  in  chemical  concentrations  over 
time  are  caused  by  the  processes  of  convective  transport,  hydrodynamic  dispersion,  and 
mixing  from  fluid  sources.  This  model  couples  the  groundwater  flow  equation  with  solute 
transport  equations. 

The  flow  equation  can  be  approximated  by  an  implicit  finite-difference  equation.  The 
model  area  is  discretized  into  a  rectangular  grid  with  each  square  being  a  node.  The  finite 
difference  equation  is  solved  numerically  for  each  node  in  the  grid  using  an  iterative 
alternating-direction  implicit  (ADI)  procedure. 

After  the  hydraulic  head  distribution  is  calculated,  the  velocity  of  groundwater  flow  can  be 
computed  at  each  node.  The  expression  for  average  velocity  of  groundwater  can  be 
derived  from  Darcy's  law.  The  groundwater  velocity  at  each  node  is  calculated  utilizing  an 
explicit  finite-difference  approximation  of  Darcy's  law.  The  computer  program  uses  an 
alternating-direction  implicit  procedure  to  solve  a  finite-difference  approximation  to  the 
groundwater  flow  equation,  and  it  uses  the  method  of  characteristics  (MOC)  to  solve  the 
solute  transport  equation.  MOC  uses  a  particle  tracking  procedure  to  represent  convective 
transport  and  a  two-step,  explicit  procedure  to  solve  a  finite-difference  equation  that 
describes  the  effects  of  hydrodynamic  dispersion,  fluid  sources  and  sinks,  and  divergence 
of  velocity. 

A  number  of  assumptions  are  inherent  in  the  solute  transport  model: 

1.  Darcy's  law  is  valid  and  hydraulic  head  gradients  are  the  only  significant  driving 
mechanism  for  fluid  flow. 

2.  The  porosity  and  hydraulic  conductivity  of  the  aquifer  are  constant  with  time,  and 
porosity  is  uniform  in  space. 

3.  Gradients  of  fluid  density,  viscosity,  and  temperature  do  not  affect  the  velocity 
distribution. 

4.  No  chemical  reactions  occur  that  affect  the  concentration  of  the  solute,  the  fluid 
properties,  or  the  aquifer  properties. 

5.  Ionic  and  molecular  diffusion  are  negligible  contributors  to  the  total  dispersive  flux. 

6.  Vertical  variations  in  head  and  concentrations  are  negligible. 
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7.  The  aquifer  is  homogeneous  and  isotropic  with  respect  to  the  coefficients  of 
longitudinal  and  transverse  dispersivity. 


Transport  Model  input 

A  model  grid  of  32  columns  by  19  rows  with  a  250  foot  lateral  spacing  was  used. 
Specified  head  cells  were  used  at  nodes  corresponding  to  the  South  Branch  of  the  Thunder 
Bay  River,  at  nodes  along  the  eastern  boundary  of  the  model  grid  area,  and  also  at  the 
sinkhole  in  the  northeastern  portion  of  the  model.  Groundwater  elevations  measured 
during  September,  1993  were  used  as  initial  input  into  the  transport  model. 

Hydraulic  conductivity  values  were  calculated  from  slug  tests  performed  at  Alpena  CRTC 
(Engineering  Science,  1989;  Earth  Technology,  1994).  Values  of  hydraulic  conductivity 
range  from  12  feet/day  at  Site  4  to  278  feet/day  at  Site  3. 

Aquifer  thickness  values  were  obtained  from  drilling  records  of  monitoring  wells  and  soil 
borings  obtained  from  the  SI  and  Rl  field  activities.  Values  listed  are  from  logs  in  which 
the  thickness  of  the  shallow  aquifer  was  clearly  discernible,  and  ranged  from  20  feet  at 
Site  5  to  65  feet  at  Site  8.  Transmissivity  values  were  calculated  by  multiplying  the 
calculated  hydraulic  conductivity  values  by  the  aquifer  thickness.  Transmissivity  ranges 
from  420  ftVday  at  TF4-MW3  to  1  5,290  ft^/day  at  CG3-MW5. 

Monitoring  of  the  discharge  of  springs  into  the  sinkhole  was  performed  during  the  SI 
(Engineering  Science,  1990)  and  an  estimate  of  approximately  18,000  gallons  of  water  per 
day  discharging  into  the  sinkhole  was  calculated.  In  order  to  obtain  a  numerical  estimate 
of  discharge  into  the  sinkhole  for  the  model,  MODFLOW  (McDonald  and  Harbaugh,  1988), 
a  3-  dimensional  finite-difference  groundwater  flow  model  was  used.  MODFLOW  was 
used  because  of  its  ability  to  simulate  the  effect  of  head-dependent  groundwater  flow  into 
a  groundwater  sink  (i.e.  the  sinkhole).  This  package  was  not  available  in  MOC.  The  same 
model  parameters  and  boundary  conditions  were  used  within  MODFLOW  as  in  MOC. 

Based  upon  hydraulic  head  data  collected  in  September  1993,  discharge  from  the  shallow 
aquifer  into  the  sinkhole  is  approximately  30,000  gallons  per  day. 

The  dispersivity  of  an  aquifer  in  two  dimensions  is  described  by  the  longitudinal  dispersion, 
the  transverse  dispersion  and  the  ratio  of  the  two  (Fetter,  1993).  As  a  contaminant  plume 
moves  further  from  its  initial  location  within  the  aquifer  by  advection  with  the  groundwater 
flow,  the  plume  spreads.  The  spreading  in  the  direction  of  groundwater  flow  is  the 
longitudinal  dispersion,  the  spreading  in  the  direction  perpendicular  to  the  groundwater 
flow  is  known  as  the  transverse  dispersion  (Fetter,  1993).  The  values  of  the  dispersion 
coefficients  are  typically  derived  via  bench  scale  tests,  aquifer  tests,  or  calibration  of 
contaminant  transport  models.  Since  no  data  presently  exists  describing  dispersivity 
within  the  shallow  aquifer  beneath  the  Alpena  CRTC  and  insufficient  data  exists  to  allow 
for  derivation  of  dispersivity  via  model  calibration,  moderate  values  of  100  feet  for 
longitudinal  dispersivity  and  30  feet  for  transverse  dispersivity  were  chosen  (Gillham  and 
Cherry,  1982).  A  more  complete  description  of  the  model  is  given  in  the  report. 

Preliminary  Groundwater  Modeling  Effort,  Earth  Technology,  August  1993. 
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Model  Calibration: 


The  groundwater  flow  model  was  calibrated  with  respect  to  the  September  1993 
groundwater  elevation  measurements.  Calibration  of  the  groundwater  flow  model  was 
accomplished  by  defining  a  set  of  parameters,  boundary  conditions,  and  stresses  that 
produce  simulated  heads  and  fluxes  that  match  field-measured  values  within  a 
preestablished  range  of  error  (Anderson  and  Woessner,  1992).  In  order  to  match  field 
measured  values  for  hydraulic  head  as  determined  during  September  1993,  a  few 
modifications  were  made  to  the  preexisting  groundwater  flow  model.  These  changes 
included  updating  the  initial  head  array,  modeling  the  sinkhole  as  a  constant  head  cell  to 
account  for  the  large  gradient  changes  in  the  vicinity  of  the  sinkhole  and  including  recharge 
to  the  model  at  a  rate  of  9  inches  per  year  over  the  whole  model  area.  By  adjusting  these 
parameters,  an  acceptable  level  of  calibration  was  achieved.  An  acceptable  level  of 
calibration  was  defined  as  a  root  mean  squared  error  (RMS)  of  less  than  2  feet.  The  RMS, 
or  the  standard  deviation  is  the  average  squared  difference  in  measured  and  simulated 
heads  and  is  given  by  the  equation: 


RMS=[l/nY^ 

i-1 


n  =  number  of  wells 
h„,  =  measured  head 
h,  =  model  simulated  head 

Certain  portions  of  the  model  may  have  values  above  the  goal  of  2  feet  while  others  fall 
much  below  this  value.  The  RMS  represents  the  average  error  present  in  the  model. 

The  following  provides  a  summary  of  the  final  calibrated  heads  for  the  flow  model. 


Well  # 

Row,  Col 

Measured 

Head  (ft) 

Simulated 

Head  (ft) 

Measured  - 
Simulated 

Head  (ft) 

MP2MW1 

10,25 

679.69 

676.68 

3.01 

MP2MW3 

10,22 

675.51 

675.46 

0.05 

MP2MW4_5  • 

11,23 

675.34 

675.70 

-0.36 

MP2MW6 

11,20 

674.86 

674.52 

0.34 

CG3MW1 

6,24 

677.38 

676.69 

0.69 

CG3MW2 

7,22 

676.29 

675.78 

0.51 

CG3MW3 

9,23 

676.50 

675.98 

0.52 

CG3MW4_5* 

8,23 

676.41 

676.08 

0.33 
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Well  # 

Row,  Col 

Measured 

Head  (ft) 

Simulated 

Head  (ft) 

Measured  • 
Simulated 

Head  (ft) 

CG3MW7 

8,20 

675.64 

674.80 

0.84 

TF4MW1 

7,9 

667.23 

668.83 

-1.60 

TF4MW2 

8,10 

659.61 

664.22 

-4.61 

TF4MW3_4» 

9,10 

658.21 

660.06 

-1.85 

SF5MW1 

12,6 

674.15 

671.34 

2.81 

SF5MW2 

14,5 

675.32 

672.67 

2.65 

SF5MW3_4* 

13,6 

676.46 

671.64 

4.82 

SF5MW6 

13,5 

674.26 

672.40 

1.86 

LF6MW1 

14,9 

672.1 

670.40 

1.70 

LF6MW2 

14,8 

672.68 

670.72 

1.96 

LF6MW3 

14,10 

671.93 

670.21 

1.72 

LF6MW4 

16,7 

672.75 

671.91 

0.84 

LF6MW5 

15,7 

673.07 

671.56 

1.51 

LF6MW6 

13,10 

671.17 

669.67 

1.50 

LF6MW8 

15,9 

673.12 

670.78 

2.34 

HN8MW1 

5,22 

676.96 

675.93 

1.03 

HN8MW2 

6,19 

675.31 

674.50 

0.81 

HN8MW3_4* 

7,21 

676.01 

675.35 

0.66 

RT9MW1 

6,16 

673.06 

672.78 

0.28 

RT9MW2 

7,14 

668.21 

670.81 

-2.60 

RT9MW3 

9,15 

670.72 

671.26 

-0.54 

RT9MW4_5* 

8,14 

667.47 

670.32 

-2.85 

RT9MW6 

8,16 

670.58 

672.33 

-1.75 

S1MW2 

13,26 

677.39 

676.63 

0.76 

S1MW3 

13,27 

677.15 

676.98 

0.17 

S1MW1  1 

15,24 

675.72 

675.67 

0.05 

S1MW12 

16,25 

674.55 

675.94 

-1.39 

I 
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Weil  # 

Row,  Col 

Measured  - 
Head  (ft) 

Simulated 

Head  (ft) 

Measured  - 
Simulated 

Head  (ft) 

S1MW13 

14,24 

675.21 

675.76 

-0.55 

S1MW14 

14,25 

673.92 

676.14 

-2.22 

MP2MW2 

12,24 

675.57 

675.96 

-0.39 

*  Indicates  that  more  than  one  well  is  present  in  each  node  and  an  average  value  for 
hydraulic  head  \was  used. 

Sum  of  Squared  Residuals  =  128.78/38  =  3.3891 
Root  Mean  Squared  Error  =  1.84 


It  should  also  be  noted  that  the  model  was  calibrated  with  respect  to  the  September  1993 
water  level  data  and  should  only  be  considered  calibrated  with  respect  to  this  data.  More 
information  on  the  water  level  fluctuation  through  time  would  be  needed  to  perform  a 
transient  calibration.  The  model  was  not  calibrated  with  respect  to  concentration  data,  but 
only  with  respect  to  hydraulic  head. 

Model  Assumptions  and  Limitations 

*  The  model  domain  consisted  only  of  the  shallow  unconfined  aquifer  (i.e.  one  layer). 

•  The  initial  head  data  input  to  the  transport  model  are  results  of  measurements  taken  in 
September  1993. 

*  Initial  concentrations  of  compounds  are  results  of  the  Round  IV  sampling  event  which 
was  conducted  from  July  to  September  1993. 

•  Hydraulic  conductivity  values  are  the  result  of  slug  tests  performed  in  November,  1987 
and  September  1993. 

*  The  model  was  calibrated  with  respect  to  hydraulic  head  using  September  1993  water 
level  data  and  should  only  be  considered  calibrated  with  respect  to  September  1993 
water  level  data. 

•  The  flow  model  was  assumed  to  be  at  steady-state  with  respect  to  hydraulic  head. 


Site  9  /  gwmod.sS 


W-22 


Sitg.  9 


Several  chemicals  were  detected  at  Site  9  above  MDNR  Type  A  or  Type  B  criteria.  These 
chemicals  include  PCE,  Benzene,  1,4-  Dichlorbenzene,  Lead,  and  2-  Methyinapthalene  at 
concentrations  of  1.5,  3.9,  18,  15.9,  and  47  ug/l,  respectively.  These  chemicals  were  all 
present  at  well  RT9MW6,  while  PCE  was  detected  In  wells  RT9MW4  and  RT9MW5. 
However,  due  to  numberical  errors  associated  with  the  very  large  gradient  changes  near 
the  sinkhole,  PCE  was  only  modeled  in  RT9MW6.  The  chemical  data  was  input  to  the 
model  at  the  node  corresponding  to  well  RT9MW6  and  concentrations  for  the  various 
analytes  were  monitored  in  the  sinkhole  with  respect  to  time. 
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